AO-A105  604  FOREIGN  TECHNOLOGY  DIV  WRIGHT “PATTERSON  AFB  OH  F/6  5/2 

JOURNAL  OF  ENGINEERING  THERMOPHYSICS  (SELECTED  ARTICLES! # (U) 

SEP  61  R  YAN#  Z  QIAN#  J  XU#  Z  JIANG#  J  YANG 
UNCLASSIFIED  FTD“lD(RS)T-0570-ei  NL 


FTD-ID(F.C)T-0570-8l 


FOREIGN  TECHNOLOGY  DIVISION 


JOURNAL  OF  ENGINEERING  THERMOPHYSICS 
(Selected  Articles) 


Approved  for  public  release 
distribution  unlimited. 


81  1  0  7  1  ^ 


FTD-ID(RS)T-0R70-81 


EDITED  TRANSLATION 


]4  y>~~FTD-^ID(RS  )T-/)57  0-8~l~) 

MICROFICHE  NR:  FTD-81-C-000846 


w - 

/  16  Sept 


{£>y  JOURNAL  OF  ENGINEERING  THERMOPHYSICS 
i  ^Selected  Articles)-  - — - - 


English  pages.: _ 7 6 


>  ■  v  ^  *  >  ■' 


j  .-'i  n 


Source: 

Country  of  origin;/  (China)  yy'. 


Fttan  1, 


Country  of  origin  j  (ChinaJ  yy.  ,»..y  y 
Translated  by:  SCxTR'AN'" 

F33657-7  8-D-0619 

Requester:  FTD/TQTA  ^ - 

Approved  for  public  release;  distribution 
unlimited.  - - -  - -  - 

f\U-  )  i  '^ 

-  /  .  /  (jj  I  \  ,  ! 

jx/'t-i V~,Ayi  >  f 

fj*  .  ■'  '  /  Xu  j 


,  V‘ 


V  ■ 

FcL  f/. 


- ---L.  - ~ - >  * 

THIS  TRANSLATION  IS  A  RENDITION  OF  THE  ORIGI- 

. i . 

NAL  FOREIGN  TEXT  WITHOUT  ANY  ANALYTICAL  OR 
EDITORIAL  COMMENT.  STATEMENTS  OR  THEORIES 

PREPARED  BY: 

ADVOCATED  OR  IMPLIED  ARE  THOSE  OF  THE  SOURCE 
ANODO  NOT  NECESSARILY  REFLECT  THE  POSITION 

TRANSLATION  DIVISION 

OR  OPINION  OF  THE  FOREIGN  TECHNOLOGY  Dt- 

FOREIGN  TECHNOLOGY  DIVISION 

VISION. 

WP.AFB,  OHIO. 

FTD-id(rs)t-or 70-81 


Date  16  Sen  19  8- 


TABLE  OF  CONTENTS 

Theoretical  Investigations  and  Experimental  Researches 

for  Higher  Subsonic  Two-Dimensional  Compressor  Cascade, 

by  Yan  Au-qun  and  Qian  Zhao-yan .  1 

A  Theoretical  Analysis  of  the  Stream  Surface  of 
Revolution  With  Supersonic  Inlet  Flow  in  a  Transonic 
Axial  Compressor,  by  Xu  Jian-zhong,  Jiang  Zuo-ren, 

Yang  Jin-sheng,  Zhang  Ying  and  Du  Zhuan-she .  35 

Preliminary  Research  on  Ablation  Experiment  Techniques 

in  Turbulent  Flow  Over  a  Flat  Plate  Specimen,  by 

Han  Yin-da,  Wang  Ke-xiang,  Yu  Guang-rong .  51 

A  New  Method  of  Cooling  Turbine  Vanes,  by  Gu  Wei-zao, 

Zhang  Yu-ming  and  Xu  Hong-kun .  65 


Research  on  the  Research  on  the  Recirculating  Flow  Rate 
in  the  Wake  of  Two-Dimensional  Flame  Stabilizers  (Under 
Cold  Conditions),  by  Zhang  Qing-fan . 


72 
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Yan  Ru-qun  and  Qian  Zbao-yan 

(Engine  Department,  North-Western  Polytechnical  University) 


SUMMARY 


The  principal  objective  of  this  article  is  to  make  use  of 
compressible  turbulence  flow  boundary  layer  theory  and  several 
characteristic  boundary  layer  parameters  used  in  the  theory  of 
boundary  layers,  such  as  momentum  thickness,  form  factors,  energy 
factors  and  other  combinational  parameters  of  a  similar  sort,  as 
well  as  making  use  of  simplified  gas  flow  speed  distributions  in 
the  boundary  layers  in  the  wake  behind  the  blade  cascade  and 
deducing,  on  the  basis  of  the  research  foundation  accumulated  in 
the  past  about  the  theory  of  blade  cascade  losses,  a  simplified 
equation  for  the  calculation  of  total  two-dimensional  blade  cascade 
pressure  losses  for  subsonic  and  transonic  speeds.  By  using  this 
equation,  it  is  possible  to  make  calculations  of  the  two-dimensional 
blade  cascade  pressure  losses  designed  operational  configurations 
and  off-designed  operational  configurations.  In  this  article,  the 
authors  also  introduce  a  new  graphical  method  for  determining  the 
corrected  gas  compressibility  effect  in  the  blade  cascade  of  com¬ 
pressors  under  operating  conditions  involving  higher  subsonic  and 
transonic  speeds. 

On  the  basis  of  the  formula  which  was  deduced  for  the  calcula¬ 
tion  of  blade  cascade  losses,  the  authors  carry  out  calculations  of 
blade  cascade  losses  for  given  blade  cascade  blade  configurations 
under  different  operational  conditions.  Besides  this,  for  given 
blade  cascade  blade  configurations,  the  authors  also  carry  out  air 
flow  tests  in  order  to  measure  the  blade  cascade  losses  under 


This  article  was  read  at  the  3rd  Annual  All-China  Thermophysics 
Engineering  Technology  Conference  in  Guilin  in  April  of  1980, 


different  conditions.  After  these  two  types  of  data  are  obtained, 
a  comparision  is  carried  out  between  them  and  this  comparison 
reveals  that  the  two  types  of  data  are  in  relatively  good  agree¬ 
ment  . 

THE  DEDUCTION  OF  THE  BASIC  EQUATION 

This  article  operates  on  the  foundation  of  research  presented 
in  references  [1,2,3]  concerning  the  theory  of  the  blade  cascade 
losses  in  low  and  high  subsonic  speed  compressors  and  it  makes  use 
of  the  theory  of  compressible  turbulent  flows  and  the  basic  prin¬ 
ciples  of  flow  through  blade  cascades,  that  is,  principles  such  as 
the  conservation  of  energy,  continuity  equations,  and  so  on,  in 
order  to  deduce  a  formula  for  the  calculation  of  the  parameter  of 
total  blade  cascade  pressure  losses  5J  for  comparision  operating  at 
subsonic  and  transonic  speeds.  In  the  process  of  deducing  this 
equation,  appropriate  use  was  made  of  several  assumptions.  The 
total  pressure  losses  for  air  flow  through  blade  cascades  is  figured 
by  making  use  of  the  average  quantity  of  total  pressure  loss. 

If  one  takes  a  look  at  Figure  1,  one  will  see  that  we  have 
determined  that  the  measurement  station  at  the  intake  of  the  blade 
cascade  is  cross-section  1-1.  The  place  in  the  flow  below  the 
blade  cascade  at  which  the  arc  lengths  are  1/2^1  is  the  measurement 
station  at  the  exhaust  of  the  blade  cascade  which  is  cross-section 
2-2.  The  average  amount  of  total  pressure  loss  which  one  finds 
between  the  cross-section  at  the  intake  of  the  blade  cascade  and 
the  cross-section  at  the  exhaust  of  the  blade  cascade  is 


In  this  equation,  P*  is  the  total  pressure  of  the  gas;  p  is  the 
density  of  the  gas;  y  is  the  tangential  direction  of  the  blade 
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cascade;  a  is  the  axial  direction  of  the  blade  cascade;  and  t  is 
the  interval  distance  of  the  blade  cascade. 


We  make  the  assumption  that  the  total  pressure  of  free  flow 
at  the  intake  of  the  blade  cascade  is  equal  to  the  total  pressure 
of  the  exhaust  of  the  cascade,  that  is  to  say,  p*  —  P7,.u  .  The 

angle  below,  fs,  represents  the  free  flow,  and  on  the  basis  of 
this  representation,  we  can  obtain 


•  t 

I 


From  Bernoulli’s  equations,  P/*.»  “  pu,i  +  +  et)  and 

P*  “*  Pi +  (l/2)’pjFi(l  +  e2)  .  in  these  equations, 

—  (1/4)  •  M2  +  1/40A/ 1  +  •  •  • ,  .  m  is  the  gas  flow  intake  Mach  number 

in  front  of  the  cascade. 


We  also  make  the  assumption  that  ph.i  pi  ,  and  on  this 
basis,  we  obtain 


p?.a  -  pi  -  i/2tpu*n*  +  tw\ 

+  Cfiiu  —  PjV  je2)  ] 


(3) 


Figure  1.  A  schematic 
diagram  of  the  boundary 
layer  distribution  along 
the  surface  of  blades  and 
of  the  pressure  distribu¬ 
tion  in  the  wake  area 
below  the  measurement 
stations  in  the  lower  part 
of  the  flow  of  the  blade 
cascade . 

1 — boundary  layer 
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Because  of  the  fact  that  and  e2  *n  the  equations  above  are 
both  positive  values  and,  in  numerical  terms,  are  both  much,  much 
smaller  than  1,  in  order  to  simplify  the  presentation  vie  will  take 
Equation  (3)  and  drop  from  consideration  the  two  quantities  in  the 
small  parentheses.  On  the  basis  of  calculations,  the  error  which 
is  caused  in  this  way  is  relatively  small.  Proceeding  in  this  way, 
we  obtain 

(in 

n.i  -  P*  -  !/2( Pu.tv)Jtl  -  Plv\) 


If  we  take  Equation  (4)  and  substitute  it  into  Equation  (2), 
and  we  then  do  a  certain  amount  of  inferring  while,  at  the  same 
time,  we  make  use  of  our  basic  knowledge  concerning  the  theory  of 
boundary  layers,  we  can  then  obtain 
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In  this  equation,  0  is  the  momentum  thickness  of  the  boundary  layer. 
6*  Is  the  displacement  thickness  of  the  boundary  layer.  6**  Is  the 
energy  thickness  of  the  boundary  layer. 


We  can  write  out  the  continuity  equations  for  the  1-1  cross- 
section  of  the  intake  of  the  blade  cascade  and  the  2-2  cross-sec¬ 
tion  of  the  flow  below  the  exhaust  of  the  blade  cascade,  that  is 

t0  v.Fjsinfi-  (T  p,F2s5n/My 

Co) 

In  this  equation  B  is  thw  flow  angle.  If  we  take  both  sides  of 
the  equation  above  and  eliminate  PiV,„i$inPi  ,  and  then  go  through 
some  simple  deduction,  we  can  arrive  at 


v„ ,  -  ( i  -  •1Y1 

Pi  sin  Pi  V  *  /* 


(7) 


If  we  take  Equation  (7)  and  substitute  it  into  Equation  (5), 
taking  into  consideration  the  fact  that  pla**  pi  >  and  we  go 
through  a  certain  amount  of  deduction,  then  we  can  obtain 


__  mar.'  --  ■  .i  l  ■  MMi  f  aT  rfc.  £. 


(A/’*),  -  A  n\\r>  \l  (jhjYh  +  (fiy/f), 

2  VjAsin^/  ii  —  (<s;//;j’ 


The  qualitative  average  total  pressure  loss  (a/*),  of  the 
blade  cascade,  as  it  is  represented  by  Equation  (8)  above  very 
obviously  shows  the  fact  that  it  is  a  function  of  the  characteris¬ 
tic  parameters  of  the  boundary  layer  wakes  in  the  flow  below  the 
measurement  stations  on  the  trailing  edge  of  the  blade  cascade 
exhaust . 


In  order  to  take  the  various  types  of  characteristic  wake 
parameters  which  are  included  in  Equation  (8),  and  by  definition, 
put  them  into  a  perpendicular  axial  plane  so  that  they  are  perpen¬ 
dicular  to  the  direction  of  the  flow,  it  is  possible  for  one,  in 
this  case,  to  take  these  various  characteristic  parameters  in  Equa¬ 
tion  (8)  and  multiply  them  by  the  sine  of  the  included  angle  between 
the  planes  in  the  direction  of  the  free  flow  and  the  axial  direction, 
and  so  obtain  the  desired  result.  At  the  same  time,  one  Inserts  In 
Equation  (8)  a  factor  reflecting  the  configuration  of  the  boundary 
layer  //  ■=•  sm/d  and  an  energy  factor  K-‘6*+/q  ,  The  reason  for 

this,  the  gas  separation  phenomenon  which  occurs  in  turbulent  flow 
boundary  layers,  and  this  effect  is  controlled  by  changes  in  the 
numerical  values  of  H  and  K.  The  conditions  for  this  type  of 
separation  are  very  close  to  M>2.0.  and  A'  1.54,4;  .  Besides 
this,  we  can  make  use  of  the  blade  cascade  density  a  —  b/t  and  the 
arc  length  b  in  order  to  express  the  Interval  distance  of  the  blade 
cascade,  t.  If  we  then  take  the  relationships  set  out  above  and 
substitute  them  Into  Equation  (8),  we  can  obtain 


P*)  _  dZi  fd\  ( s*n  P\  V  (  JL\ 

2  Wiring  )\bh 


1  +  K, 


l  '  b  h  *ln  fit  J 


Concerning  the  comparison  of  values  which  this  article  makes 
between  the  blade  cascade  total  pressure  losses  (AP#^  and  the  blade 
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1. 

'2 


cascade  intake  dynamic  pressureAp^V^  as  measured  at  a  station 
below  the  cascade,  this  comparison  can  be  used  as  the  total  blade 
cascade  pressure  loss  coefficient  co  in  order  to  represent  the  flow 
losses  in  the  blade  cascade.  At  the  same  time,  the  following 
definitions  pertain.  02  is  the  momentum  thickness  parameter  of 
boundary  layers,  that  is  to  say,  Oi  —  (0/£)2  •  a!  sin  ft. 

If  one  takes  the  two  concepts  mentioned  above  and  substitutes  them 
into  Equation  (9),  then  it  is  possible  to  obtain  total  pressure 
loss  coefficient  of  the  blade  cascade  involved,  which  is 


a 


Ex  ( sln  ft  V  § .  1 

Pj  'sin  ft/  2  ll—dj-ItY 


(10) 


In  order  to  calculate  the  value  of  u  in  Equation  (10),  it  is 
necessary  to  establish  governing  principles  for  a  flow  model  of  the 
flow  speed  distribution  in  the  wake  of  the  measuring  stations  below 
the  exhaust  of  blade  cascades.  In  order  to  accomplish  this,  we 
make  use  of  a  simplified  form  of  gas  flow  speed  distribution  model 
for  wakes,  as  shown  in  Figure  2.  In  Figure  2,  the  vertical  coordi¬ 
nates  represent  the  speed  ratio  V/Vfs,  and  the  horizontal  coordi¬ 
nates  represent  the  boundary  layer  thickness  ratio  y/(stmlt/ 2).  In 

this  equation,  5  is  the  boundary  layer  thickness  on  the  surface  of 
the  blades,  and  8i»n  is  the  total  boundary  layer  thickness  above 
and  below  the  blade  surfaces.  As  far  as  the  simplified  gas  flow 
speed  distribution  is  concerned,  it  is  possible  to  write  out  the 
form  of  the  relationship  below,  that  is 

v/vu  -  <y/(  W 2))* 


Figure  2.  Simplified 
gas  flow  speed  distribution 
model  for  wake. 

1 — speed  ratio 
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In  this  equation,  n  is  an  exponent.  If  we  are  considering  the 
case  of  the  blade  cascades  of  compressors  in  subsonic  and  tran¬ 
sonic  applications,  then  n  =  0.2^0. 4.  Besides  this,  if  we  are 
speaking  of  the  relationship  between  the  energy  factor  K  and  the 
configuration  factor  H  as  they  appear  in  this  type  of  simplified 
flow  speed  distribution,  and  if  we  consult  references  [1,2,3]  in 
order  to  make  use  of  what  they  present,  then  the  relationship  is 
K  —  (H  +  l)/(3f/  —  1)  .  If  we  then  take  this  relationship  and 

substitute  it  into  Equation  (10),  it  is  possible  to  obtain 


o  -  (iisAV 02( i  _  QH  )-i  .,..1%...  . 
•  PjVsln  fij  3H,-I 


(ID 


What  one  sees  above  is  nothing  more  than  a  formula  for  calcul¬ 
ating  the  blade  cascade  loss  coefficient,  DJ,  for  subsonic  and  tran¬ 
sonic  compressors  as  that  formula  was  deduced  from  the  theory  of 
compressible  turbulent  flow  boundary  layers.  It  can  be  seen  from 
Equation  (11)  that  this  method  of  calculation  can  clearly  be  used 
for  carrying  out  calculations  of  total  pressure  loss  coefficients 
in  blade  cascades  when  those  cascades  are  operating  in  a  slow  speed 
operational  configuration,  such  as  (M^  <  0,3).  The  reason  for  this 
Is  the  fact  that,  under  such  conditions,  it  Is  possible  to  recog¬ 
nize  the  gases  Involved  as  being  non-compressible ,  that  is  to  say, 
p-^  =  c>2'  this  way.  Equation  (11)  becomes 


Actually,  Equation  (12)  Is  nothing  else  than  the  theoretical 
relationship  between  the  total  blade  cascade  pressure  losses  for  low 
rotation  speed  compressors  as  deduced  on  the  basis  of  a  non-com¬ 
pressible  flow  by  a  consideration  of  Lieblein,  S.,  and  others  work¬ 
ing  in  a  similar  field  as  found  in  [2],  Besides  this.  It  can  be 
pointed  out  on  the  basis  of  references  [5,6]  that  when  one  Is  deal¬ 
ing  with  a  situation  in  which  a  blade  cascade  is  at  a  high  subsonic 
speed  so  that  the  Mach  number  even  up  to  the  intake  in  front  of  the 
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cascade  reaches  1  or  1.2,  then  the  blade  trough  shock  loss  coeffi¬ 
cient  od  (no  matter  whether  we  are  talking  about  blunt  blades  or 
sharply  pointed  blades)  is  relatively  small  and,  in  general,  does 
not  exceed  an  approximate  value  of  0.01.  Because  of  this,  gener¬ 
ally  speaking,  when  one  gets  involved  with  the  calculation  of 
blade  cascade  losses  in  situations  in  which  the  rotation  speeds 
involved  are  high  subsonic  to  transonic,  it  is  possible  to  almost 
leave  out  shock  losses  inside  the  blade  cascade  from  any  consider¬ 
ation.  Just  as  is  pointed  out  in  [7],  when  the  intake  angle  of  the 
blade  cascade  is  increased,  the  blade  cascade  losses  must  also  be 
increased;  however,  at  such  a  time,  the  main  reason  for  this  pheno¬ 
menon  is  not  an  increase  in  shock  losses  in  the  blade  paths,  but 
rather  is  primarily  due  to  an  increase  in  the  blade  form  losses 
The  reason  for  this  is  the  fact  that  the  magnitude  of  shock  wave 
losses  in  blade  cascade  paths,  for  all  ranges  of  change  in  intake 
angle,  do  not  exhibit  obvious  changes.  Of  course,  it  is  also 
possible,  on  the  basis  of  the  method  discussed  in  [8]  to  make  a 

simple  calculation  of  the  coefficient  of  shock  wave  losses  to  under 

s 

transonic  operational  conditions  in  the  blade  cascade. 

Overall,  what  was  said  above  is  this:  Concerning  Equation  (11) 
which  was  deduced  in  this  article,  it  is  possible  to  use  it  in  the 
calculation  of  the  total  pressure  loss  coefficients  of  compressor 
blade  cascades  (whether  one  is  dealing  with  blunt  blades  or  sharply 
pointed  blades)  under  designed  and  non-designed  operational  condi¬ 
tions  of  low  speeds,  high  subsonic  speeds  and  transonic  speeds.  The 
reason  for  this  is  that,  under  different  operational  conditions, 
the  magnitude  of  blade  cascade  losses  is  generally  a  reflection  of 
changes  in  the  numerical  values  of  various  types  of  boundary  layer 
characteristic  parameters  in  the  wake  of  the  cascade  involved. 

Below  we  will  discuss  the  quantity  P^/P2  as  found  in  Equation 
(11).  This  quantity  is  principally  concerned  with  a  consideration 
of  the  effects  of  gas  compressibility  on  blade  cascades  under  high 
subsonic  and  transonic  intake  speed  conditions.  If  we  give 
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consideration  to  the  relationship  between  gas  form  equations  and 
the  gas  kinetics  total  temperature  and  static  temperature  as  well 


as  total  pressure  and  static  pressure  and,  at  the  same  time, 

assume  that  P*  03  P*.i  and,  if  we  also  make  use  of  the  relation- 
£ 

ship  ,  then  on  the  basis  of  all  this,  we  can  obtain 


£\ 

Pi 


('  + 


k~  1 

2 

K-  1 
2 


* 


1 


K  - 1 


m! 


1  4-  Mi 

2 


Figure  3-  The  curve  defined 
by  the  relationship  between 
the  blade  cascade  intake  Mach 
number  coefficient  M*  and  the 
Mach  number  Mi  of  the  flow  com¬ 
ing  into  the  blade  cascade. 

1 —  sharply  pointed  blade  tip  forms; 

2 —  blunt  blade  tip  forms 

3 —  Northwestern  Folytechnical 
University  blade  cascade  data — 
Pointed  blade  form 
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If  we  take  the  various  quantities  on  the  right  side  of  the 
equation  above  and  make  a  series  expansion  of  them  while,  at  the 
same  time,  eliminating  the  M  number  quantity  for  the  high  harmonic, 
then  we  obtain 

Pi/ft  -  0  -  1/2*  Mj)(  1  +  1/2  •  Mj)  -  M,  (_i3) 

The  quantity  A?i  in  Equation  (.13)  above  is  called  the  intake 
Mach  number  coefficient.  The  numberical  value  for  the  magnitude 
of  this  quantity  is  primarily  determined  by  the  Mach  number  M^  for 
the  flow  coming  into  the  cascade  concerned.  Because  of  this,  when 
one  is  considering  the  time  period  in  which  the  high  speed  gas 
flow  is  going  through  the  blade  cascade  in  Equation  (11),  he  should 
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be  concerned  mostly  with  the  influence  of  the  compressibility  of 
the  gases  involved.  In  this  way,  under  operating  conditions 
involving  subsonic  and  transonic  speeds,  the  blade  cascade  total 
pressure  loss  coefficient  w,  can  be  written  as 

-  W)-‘*  4//«/<3/f,-  O  (1*0 

Concerning  the  blade  cascade  intake  Mach  number  coefficient 
.  A?i  as  it  appears  in  Equation  (14)  above,  it  is  possible  to  deter¬ 
mine  its  value  by  experimentation.  Because  of  this  fact,  we  took 
set  7  of  the  two  dimensional  experimental  data  for  the  blade  cas¬ 
cades  involved,  including  within  it  the  third  set  of  data  from 
[9]  (including  data  on  both  the  blunt  blade  form  and  the  pointed 
blade  form)  as  well  as  the  set  four  experimental  data  produced  by 
Northwestern  Polytechnical  University  for  the  blade  cascades 
involved  (pointed  blade  form)  [10],  and  we  used  them  in  Equation 
(13)  in  order  to  carry  out-  calculations.  In  this  way,  it  is  poss¬ 
ible  to  obtain  values  for  Mi  for  different  blade  cascade  intake 
Mach  numbers  M-^.  By  doing  things  this  way,  we  use  SJi  values  as 
!  vertical  coordinates  and  take  intake  Mach  numbers  M^  as  horizontal 

coordinates  in  order  to  draw  out  the  curve  defined  by  the  relation- 
ship  M,  —  m,  ,  as  it  is  shown  in  Figure  3.  From  Figure  3,  it  can 
be  seen  that  values  of  for  different  values  of  M^  all  fall 
within  a  relatively  narrow  band  of  data.  Because  of  the  experimental 
data  which  was  employed  to  make  Figure  3,  obtained  as  it  was  from 
flow  parameters  for  different  blade  cascade  blade  form  geometry 
parameters  and  operational  conditions,  the  curve  M,— A/,  which 
is  presented  in  Figure  3  is  of  a  type  which  can  be  used  relatively 
widely.  Concerning  the  graph  of  A?i—  Mt  which  we  have  spoken  of 
above,  it  was  made  using  a  new  type  of  method  presented  by  the  author 
and  which  corrects  for  the  effects  of  compressibility  of  gases  and 
various  types  of  operating  conditions  involving  compressor  blade  cas¬ 
cade  intake  speeds  in  the  high  subsonic  and  transonic  speed  ranges. 
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A  COMPARISON  WITH  EXPERIMENTAL  DATA 

To  start  with,  we  made  use  of  Equation  (14)  which  expresses 
a  formula  for  the  total  blade  cascade  pressure  loss  coefficient  u 
and,  on  the  basis  of  this  formula,  we  carried  out  calculations 
for  compressor  blade  cascades  using  standard  blade  forms  under 
different  types  of  operational  conditions.  The  exponent  for  the 
simplified  flow  form  model  of  gas  speed  sitribution  in  the  wake  of 
the  blade  cascade  is  n  =  0.35. 

Besides  this,  we  also  carried  out  wind  flow  experiments  on 
blade  cascades  of  the  type  using  the  standard  pointed  blade  form 
and  blade  cascades  using  the  blunt  blade  form,  under  both  designed 
and  off-design  operational  configurations.  The  wind  tunnel  exper¬ 
iments  on  blade  cascades  were  carried  out  at  Northwestern  Polytech- 
nical  University  using  its  two  dimensional  blade  cascade  wind  tunnel. 
Concerning  the  geometrical  parameters  and  construction  of  this  wind 
tunnel,  these  have  already  been  discussed  in  detail  in  [10,11]. 

The  blunt  blade  series  which  we  used  was  the  Sovie^°<!>A^0/27P45 

blade  form.  During  these  experiments,  the  range  of  blade  cascade 
intake  Mach  numbers  used  was  M1  =  0.40M3.95,  and  the  blade  cascade 
Intake  blade  arc  Reynolds  number  was  R,>  4.9  x  10*  .  Concerning 

the  experimentation  technology  of  this  wind  tunnel  in  terms  of  its 
automatic  measurement  and  automatic  recording  systems,  one  can  get 
a  detailed  view  of  all  of  them  by  consulting  [12]. 


In  Figure  4,  one  sees  presented,  for  blunt  blade  form  blade 
cascades  under  design  operating  conditions,  the  values  of  the  total 
pressure  loss  coefficient  (w caj_cuiated  ) »  wh:ich  are  obtained  by  carry¬ 
ing  out  calculations  using  Equation  (14)  and  the  values  for  the 
total  loss  coefficients  which  were  obtained  by  carrying  out  experi¬ 
mental  measurements  as  well  as  a  comparison  of  the  two  which  shows 
that  both  sets  of  values  are  In  very  good  agreement.  In  Table  2, 
we  see  presented  the  same  two  types  of  data,  that  is  (“caicuiated  ^ 
and  (^asured)  for  Po:inted  blade  blade  cascades  under  off-design 
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;ey  : 


A  comparison  of  the  values  of  u  ,  ,  .  ,  and  oj 

*  calculated  measured 

(design  configuration  of  pointed  blade  from  blade  cascade) 


i  -  +  5* 

M, 

0.4211 

o.ouyj 

0.7102 

0.8207 

0.8891 

L  =  2.05 

*  “It* 

U.07M) 

0.1002 

0.1120 

0.1522 

0.2180 

t 

0.0770 

0. 1014 

0.1084 

0.1529 

0.2110 

ABLE  2 


1 — Mach  number  of  incoming  flow;  2 — calculated;  3 — measured 
A  comparison  of  the  values  of  to  and 


(off-design  configuration  of  pointed  blade  form  blade 
cascade ) 


«-0* 

A  =  2.05 

t 

Mi 

1  “tt* 

-3  “*» 

0.2920 

0.0519 

0.0S12 

0.6083 

0.0780 

0.0772 

0.7170 

0.0881 

0.0834 

0.7720 

0.1915 

0.2380 

.  =  +7.5* 

0.299S 

0.6083 

0.7170 

0.7720 

0.8709 

a 

“it* 

0.0766 

0.1070 

0.1210 

0.2150 

0.220 

±  =  2.05 

s 

3 

“** 

0.0810 

0.1340 

0.1266 

0.2140 

0.2360 

$  =  +10° 

terns i v, 

0.4140 

0.5074 

0.6104 

0.7075 

0.8185 

A  _ 

a 

“it* 

O.U53 

0.1425 

0.1610 

0.1960 

0.2050 

±  **  2.05 

$ 

3 

“■* 

0.1240 

0.1410 

0.1907 

0.2750 

0.2596 

i  =  -2* 

'fcflESilfR  M, 

0.3930 

0.5981 

0.7054 

0.8000 

A  _ 

2- 

“it* 

0-0378 

0.0229 

0.0515 

0.0563 

2.05 

t 

3 

“at 

0.0457 

0.0595 

0.1885 

0.2280 

im  +5° 

0.4140 

0.6033 

0.7096 

0.8163 

0.9189 

A  .  - 

T 

“n* 

0.0515 

0.0592 

0.0740 

0.0885 

0.1210 

±  -  1.5 

1 

3_ 

0.0483 

0.0596 

0.0773 

0.0892 

0.1680 

Key:  1 — Mach  number  of  incoming  flow;  2 — calculated;  3 — measured 


operational  configurations  and  again,  the  two  are  in  relatively 
close  agreement.  In  Table  3,  we  see  presented  a  comparison  of  the 
same  two  types  of  data,  C^calculated}  and  for  blunt  blade 

form  blade  cascades  under  design  operational  configurations. 


Prom  the  various  types  of  data  presented  in  Tables  1  to  3,  it 
is  possible  to  see  that,  when  one  is  considering  the  type  of  situa¬ 
tion  in  which  a  blade  cascade  is  operating  within  design  parameter 
ranges,  the  calculations  which  are  made  on  the  basis  of  Equation  (1*0 


and  the  blade  cascade  total  pressure  loss  coefficients  w  which 
are  obtained  by  experimental  measurements  are  both  in  relatively 
good  agreement.  In  the  case  in  which  the  blade  cascade  under  con¬ 
sideration  is  operating  in  an  off-design  configuration  (Table  2), 
the  two  types  of  data,  <^calculated  and  ^measured ,  are  also  relatively 
close.  If  the  operational  configuration  of  the  blade  cascade  being 
considered  is  operating  too  far  off  the  design  specifications  (for 
example,  the  angle  of  incidence  is  i  =  +10°),  then  the  results 
from  using  Equation  (.1*0  to  do  calculations  of  w  lose  their  uni¬ 
formity  with  the  measured  values  for  the  same  quantity  and  there  is 
an  obvious  discrepancy  between  the  two.  The  reason  for  this  is 
that,  at  such  a  time,  the  air  flow  in  the  passages  of  the  blade  cas¬ 
cade  produce  severe  separation  and  in  the  entire  blade  trough  one 
sees  almost  complete  separation  of  gases.  At  the  same  time,  the 
boundary  characteristic  parameters  which  are  obtained  by  calcula¬ 
tions,  such  as  the  form  factor  H2  value,  already  approach  the  num¬ 
erical  magnitude  of  2. 

TABLE  3.  A  comparison  of  the  values  of  <^cadcudated  and 

^measured  (t3lunt  blade  form  blade  cascade  in  a 
given  operation  configuration). 

(The  row  of  data  marked  with  an  asterisk  is  a  collection 
of  the  wind  tunnel  test  characteristics  for  the  Soviet 
A-40  blade). 


0.3549 

0.5971 

0.7149 

0.8075 

z 

0.0194 

0.0238 

0.0306 

0.1135 

0.0241 

0.0267 

0.0297 

0.1110 

w  •£>#>»* 

0.0200 

0.0200 

0.0220 

0.0390 

Key:  1 — Mach  number  for  the  incoming  flow;  2 — calculated; 

3 — measured;  4 — collected  characteristics 

On  the  basis  of  what  has  been  discussed  above,  if  one  only 
needs  to  know  the  total  boundary  layer  thickness  at  the  mea¬ 

surement  station  down  the  flow  from  the  exhaust  of  the  blade  cas¬ 
cade,  then  it  Is  possible,  on  the  basis  of  Equation  (14)  to  calcu¬ 
late  the  total  blade  cascade  pressure  loss  coefficient  uT  of  blade 
cascades  under  any  certain  operational  configuration. 
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It  is  worth  'while  to  point  out  that  if  we  are  carrying  out 
new  design  calculations  for  compressor  blade  cascades  (this  does 
not  consider  the  question  of  whether  the  blades  involved  are  of 
the  blunt  type  or  the  pointed  type),  then  it  is  also  possible  to 
make  use  of  Equation  (14)  in  order  to  calculate  the  total  blade 
cascade  pressure  loss  coefficient  u  of  a  blade  cascade  in  any  cer¬ 
tain  operational  configuration  and  this  can  be  used  as  a  reference 
value  for  use  during  redesigning.  At  such  a  time,  it  is  only 
necessary,  on  the  basis  of  the  statistical  data  or  by  referring 
to  the  data  calculated  for  a  certain  blade  cascade  in  a  given 
operational  configuration,  to  calculate  the  total  boundary  layer 
thickness  for  the  upper  and  lower  blade  surfaces,  that  is, 
and  that  will  take  care  of  everything.  In  Figure  4  in  this  article, 
we  see  graphed  out  the  comparative  values  of  the  parameters  for  a 
blade  cascade  employing  a  certain  type  of  geometrical  parameters 
of  blade  form,  that  is,  for  both  the  pointed  blade  type  and  the 
blunt  blade  type,  while  employing  different  values  for  the  angle 
of  Incidence.  The  values  compared  are  the  total  boundary  layer 
thickness  for  the  blades  Involved,  that  is,  and  the  Interval 

distance  between  the  blades  of  the  cascade,  t.  The  ratio  of  these 
two  quantities  is  then  plotted  against  the  experimentally  measured 
values  for  the  intake  Mach  numbers,  The  graphic  representa¬ 

tion  of  this  relationship  Is  then  used  as  a  reference  In  the  anal¬ 
ysis  of  blade  pressure  losses  during  the  designing  of  blade  cascades 
for  compressors. 


Figure  4.  The  curve  defined 
by  the  relationship  between 
the  value  of  the  ratio  'Sfwn/t 
which  represents  the  total 
thickness  of  the  boundary 
layers  involved  and  the  inter¬ 
val  distance  between  the  blades 
for  varying  angles  of  intake 
incidence 

Key:  1 — V  cross-section  of  an 

operating  wheel  in  the  second 
stage  of  a  certain  compressor 
employing  pointed  type  blades  x's*** 

2 — blunt  type  blades  */»-*•  m  *-55.j «mm 

aSwWK"  5W40/27.6n« 

hfl  «■  1.3  *  “  40mm 
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CONCLUSIONS 


We  will  now  present  the  main  conclusions  which  can  be  drawn 
from  the  material  presented  In  this  article. 

1.  If  one  makes  use  of  the  theory  for  boundary  layers  in 

compressible  turbulence  flow,  employs  the  performance  parameters 
which  mainly  characterize  the  boundary  layers  in  the  wake  of  a 
blade  cascade  and  also  make  use  of  the  M  number  for  the  flow  coming 
into  the  cascade  involved  (that  is,  the  value)  as  well  as 

parameters  for  the  geometry  of  the  blade  cascade  blade  being  con¬ 
sidered,  then  it  is  possible  to  relate  these  quantities  to  the 
qualitative  average  pressure  loss  coefficient  for  the  blade  cas¬ 
cades  concerned  when  they  are  employed  In  subsonic  and  transonic 
modes . 

2.  This  article  makes  use  of  the  basic  equation  for  pressure 
losses  In  blade  cascades  in  order  to  deduce  other  theoretical  rela¬ 
tionships  which  are  used  together  with  a  large  amount  of  wind 
tunnel  test  data  on  blade  cascades  to  present  a  new  method  of  graph¬ 
ically  representing  Mx—\tx  ,  which  is  effective  for  blade  cascades 
in  compressors  operating  at  high  subsonic  and  transonic  speeds  and 
corrects  for  the  effects  of  the  compressibility  of  the  gases 
involved. 


3.  The  magnitude  of  the  total  pressure  losses  in  a  blade  cas¬ 
cade  is  primarily  dependent  on  the  intake  Mach  number  parameter 
'Mi-  (which  is  dependent  on  the  number  of  the  flow  coming  into 

the  blade  cascade)  and  on  the  value  of  the  ratio  02/b  between 
boundary  layer  momentum  thickness  and  blade  arc.  As  far  as  the 
boundary  layer  energy  factor  K2  or  the  configuration  factor  H2  go, 
they  are  secondary  in  their  influence  on  the  total  pressure  losses 
in  a  blade  cascade. 
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4.  If  we  consider  a  certain  given  blade  cascade  operating 
under  a  certain  type  of  configuration,  then  if  we  already  know 
c*>n  estimate  the  total  boundary  layer  thickness  on  the  surface 
of  blades  down  the  air  stream  from  the  exhaust  of  the  cascade  being 
considered,  then  it  is  possible  to  calculate  the  qualitative  aver- 
-e  total  pressure  coefficient  for  the  cascade  being  considered, 
that  is,  u  (including  both  design  operational  configurations  for 
the  blade  cascade  being  considered  and  off-design  operational  con¬ 
figurations  ) . 


REFERENCES 

11]  Nti'wart,  \V.  L. :  Analysis  of  Two-Dimensional  ('••ni|nessilile-Fli>w  loss  Cliuractoristies  Downstream 
of  Turbumuchine  lilude  Hows  in  terms  of  Basie  iioumlary-luyer  Characteristics,  A  AC  A  TN  3515, 

(1955). 

(2)  Lieblein,  8.  and  Roudvbu.sh.  W.  11.;  Theoretical  loss  Relation  for  low  Speed  Two  Dimensional 
Cascade  Flows,  .Y.4C.4  r.V  .ICtiJ,  (1950). 

[  3  ]  Lieblein,  8.  and  Roudehush.  W.  U.:  Low-Speed  Wake  Characteristics  of  Two-Dimensional  Cascade 
and  Isolated  Airfoil  Sections,  NACA,  TX  3771,  (1956). 

[4  ]  Bchlichting,  H. :  Application  of  Boundary-Layer  Theory  in  Turbomuchincry,  Train.  ASMS  81D, 

(1959). 

[51  Koch,  C.  C.  and  Smith  Jr..  L.  H. :  Loss  Sources  and  Magnitude  in  Axial  Flow  Compressors,  Tram. 

ASMS  Series  .1,  96.  3.  July,  (1976). 

[6]  Heilman,  H.,  Starkcn,  H.  and  Wegcr,  H.:  Ca.-eadc  Wind  Tunnel  Tests  on  Blade  Design  for 
Transonic  and  Supersonic  Compressors.  AD  6J7837,  Sept.  (1908). 

[  7  ]  Dowis,  O.  W.  and  Uartmann,  M.  J. :  Shock  Losses  in  Transonic  Compressor  Blade  Rows,  Trass. 

ASME  S3A,  (1901). 

[8]  Swan,  W.  C.:  A  Practical  Method  of  Predicting  Transonic  Compressor  Performance,  Trans. 

ASME  Series  A,  83.  1901. 

[  9  ]  Todd,  K.  W. :  An  Experimental  Study  of  Three-Dimensional  High  Speed  Air  Conditions  in  a 
Cascade  of  Axial  Compressor  Cascade,  ARC  RM  -’794,  (1954). 

[10]  Blade  Cascade  Experimentation  Team:  "Experimental  Research 
on  Plane  Blade  Cascades  in  Order  to  Determine  the  Effects  of 
Backward  Sweep",  Northwestern  Polytechnical  University, 
Scientific  and  Technological  Materials  Syllabus,  No,  63**,  1968. 

[11]  Blade  Cascade  Experimentation  Team:  "Plane  Blade  Cascade  Wind 
Tunnels",  Northwestern  Polytechnical  University  Press,  1977. 

[12]  Qian  Zhao-yen :  "Automatic  Recording  and  Data  Processing  in 
Wind  Tunnel  Test  Technology  Relating  to  Plane  Blade  Cascades", 
Northwestern  Polytechnical  University,  Scientific  and  Techno¬ 
logical  Materials  Syllabus,  No.  1076,  1979. 


16 


THEORETICAL  INVESTIGATIONS  AND  EXPERIMENTAL 
RESEARCHES  FOR  HIGHER  SUBSONIC  TWO- 
DIMENSIONAL  COMPRESSOR  CASCADE 

Yan  Ru-qun,  Qian  Zhao-yan 
(Engine  Department,  North-Western  Polytechnical  University) 

Abstract 

The  losses  occured  in  air-flow  of  compressor  cascade  are  mainly  to  be  controlled  by 
occufing  developing  and  separting  of  boundary-layer  on  the  surfaces  of  blade.  There¬ 
fore,  by  means  of  adapting  the  theory  of  compressible  turbulent  boundarylayer,  utiliz¬ 
ing  some  characteristic  parameters  such  as  momentum  thickness,  form  factor,  energy 
factor  and  other  combination  parameters,  and  then  supposing  the  flow  model  of  simple 
power  of  velocity  distribution  in  the  wake  behind  the  cascade,  a  simplified  calculating 
equation  for  total  pressure  loss  coefficient  of  subsonic  and  transonic  cascades  of  the 
compressor  is  derived.  The  basic  calculating  equation  derived  would  be  used  for  predic¬ 
ting  the  total  pressure  losses  of  both  design  condition  and  off-design  condition!  of  the 
compressor  cascade. 

In  this  paper,  a  new  curvilinear  method  for  correcting  the  gaseous  compressible 
effect  under  higher  subsonic  and  transonic  intake  condition  upstream  of  the  cascade  is 
presented. 

According  to  this  basic  equation,  the  total  pressure  losses  for  given  compressor 
cascade  are  calculated,  and  then  tested  on  the  plane  wind  tunnel  at  various  conditions, 
both  the  results  are  in  better  agreement. 
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THE  SHOCK  WAVE  STRUCTURE  AUD  THE  ANALYSIS  OF  THE  FLOW  THROUGH 
TRANSONIC  TURBINE  CASCADES* 


Min  Da-fu 

(The  Qing  Ho  Aerodynamic  Lab.  of  the  Beijing  Heavy  Electrical 
Machine  Works) 

SUMMARY 

On  the  basis  of  Schlieren  photographs  which  were  taken  during 
tests  run  on  plane  blade  cascade  air  flow  during  transonic  turbine 
operation,  an  analytical  discussion  is  carried  out  on  the  subject 
of  shock  wave  structure  and  development  in  exhaust  flow  fields  and 
blade  cascade  paths.  On  the  basis  of  the  results  of  these  experi¬ 
ments,  we  verified  the  existence  of  shock  wave  deflection  on  free 
boundary  surfaces  as  well  as  the  influence  of  this  phenomenon  on 
the  convection  flow  characteristics.  In  the  third  section  of  this 
article,  we  carry  out  a  discussion  of  several  classic  types  of  the 
phenomenon  of  interference  between  shock  waves  and  boundary  layers. 

1 .  FOREWORD 

The  shock  wave  system  structure  in  the  passages  of  plane  blade 
cascades  as  used  in  transonic  turbine  applications  and  in  exhaust 
flow  fields  as  well  as  the  interference  between  these  shock  waves 
and  the  boundary  layers  involved  determine,  to  a  very  large  extent, 
the  flow  characteristics  of  a  blade  cascade.  In  order  to  do  research 
on  the  design  principles  of  transonic  blade  cascades  for  the  purpose 
of  developing  blade  cascades  which  do  not  produce  shock  waves,  it  is 
necessary  to  make  very  careful  observations  and  analyses  of  the 
shock  wave  structures  involved  in  winding  flow  fields.  Our  labor¬ 
atory  carried  out  experiments  on  several  sets  of  blade  cascades  as 


*  This  paper  was  read  at  the  Third  Annual  All-China  Engineering 
Thermophysics  Technology  Conference  in  Gui  Zhou,  April  1980. 
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well  as  taking  Schlieren  photographs  to  correspond  with  these  cas¬ 
cades  during  the  course  of  the  experiments.  The  results  of  the 
experiments  showed  that  when  the  exhaust  M  number  approached  1, 
that  shock  waves  appeared  in  the  passages  of  the  blade  cascades 
involved  as  well  as  in  the  exhaust  flow  fields.  When  this  M  num¬ 
ber  was  extended  somewhat  more,  the  strength  of  the  shock  wave 
involved  as  well  as  their  positions  exhibited  changes.  Reference 
[1]  has  already  documented  and  discussed  these  problems  and  this 
article  will  take  the  discussion  of  models  for  the  shock  wave  systems 
under  consideration  a  bit  further  on  the  basis  of  the  results  of 
experimentation. 

Concerning  plane  blade  cascade  wind  tunnels,  due  to  the  fact 
that  the  blades  are  ve.-’y  limited,  in  the  flow  from  the  first  blade 
to  the  last  one,  there  exist  on  the  boundaries  of  the  flow  free 
boundary  surfaces.  The  outside  tail  shock  waves  produced  by  the 
trailing  edges  of  the  blade  cascades  being  considered  will  produce, 
on  the  free  boundaries,  deflections  and  the  formation  of  complex 
systems  of  shock  waves.  Varying  with  changes  in  the  M  number, 
these  types  of  complex  structures  will  also  exhibit  changes.  [2J 
has  already  given  its  main  attention  to  this  question.  It  can  be 
verified  from  the  Schlieren  photographs  of  the  several  blade  cascades 
which  went  through  testing  in  our  laboratory  that  this  sort  of 
phenomenon  does  exist  and  we  will  set  out  from  the  basic  principles 
of  gas  flow  dynamics  to  arrive  at  a  solution  of  this  problem  area. 

The  essence  of  this  problem  is  that,  due  to  the  deflection  which 
takes  place  with  this  type  of  shock  wave,  one  gets  the  formation  of 
multi-peak  type  parameters  and  wave  forms  in  the  exhaust  flow  fields 
of  turbine  blade  cascades  in  transonic  applications.  These  types  of 
multi-peaked  forms  also  exert  an  influence  on  our  evaluation  of  the 
aerodynamic  performance  of  blade  cascades.  This  article  will  take 
the  results  of  aerodynamic  experiments  and  the  available  Schlieren 
photography  in  order  to  analyze  and  solve  this  problem. 
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The  interference  between  shock  waves  and  boundary  layers  is 
an  important  factor  affecting  turbine  blade  cascade  pressure  losses 
in  transonic  applications  as  well  as  other  performance  character¬ 
istics.  From  Schlieren  photographs,  we  can  obtain  some  enlighten¬ 
ment  on  the  subject  of  how  to  solve  the  problem  of  this  physical 
phenomenon.  However,  due  to  the  excessive  complexity  of  the  prob¬ 
lem,  it  is  difficult  even  from  the  photography  which  we  already 
have  to  reach  firm  conclusions  concerning  this  problem  and  we  still 
must  wait  upon  continued  efforts  in  the  future. 

2.  EXPERIMENTAL  EQUIPMENT  AND  METHODS 

All  the  experimentation  which  we  did  was  carried  out  on  the 
F01 intermittent  type  transonic  turbine  blade  cascade  wind  tunnel 
[31.  The  Schlieren  photography  was  done  by  making  use  of  a  505-1 
type  Schlieren  apparatus  [4]  and  the  camera  employed  was  a  type  DF-7 
Seagull  brand.  The  speed  selected  for  the  photography  was  1/500 
sec*\  During  the  time  when  the  photography  was  being  done,  the 
semi-circular  dish  in  the  experimental  section  of  the  wind  tunnel 
was  modified  to  include  organic  glass  observation  windows.  Given 
this  modification,  the  installation  plate  for  the  blade  cascade, 
due  to  strength  and  equipment  capabilities,  could  not  help  but  make 
use  of  organic  glass  in  its  manufacture  and  due  to  flaws  in  the 
optical  capabilities  of  these  windows,  the  clarity  and  precision  of 
the  images  of  the  gas  flow  in  the  passageways  of  the  blade  cascade 
involved  fall  somewhat  short  of  what  could  otherwise  be  expected. 

Based  on  the  principles  governing  the  formation  of  Schlieren 
images,  it  is  only  when  one  has  had  a  density  gradient  In  the  dir¬ 
ection  of  the  normal  line  of  the  cutting  edge  of  a  blade  that  It 
is  possible  to  get  the  production  of  a  clear  and  precise  image. 
Because  of  this,  if  we  take  the  knife  edge  and  position  it  in  a 
direction  perpendicular  to  the  gas  flow  in  the  exhaust,  then  it  will 
be  possible  to  take  more  completely  satisfactory  pictures  of  the 
shock  wave  systems  under  study.  Moreover,  when  the  knife  edge  is 
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placed  In  a  direction  parallel  to  the  direction  of  the  can  flow, 
then  it  is  possible  to  achieve  even  clearer  and  more  precise  pic¬ 
tures  of  the  wakes  present. 


3.  EXPERIMENTAL  RESULTS  AND  ANALYSIS 

1 .  The  formation  and  development  of  the  shock  wave  systems 
In  the  blade  cascade  passageways 

Figure  1  consists  of  Schlieren  photographs  of  certain  blade 
cascades  under  various  operating  conditions.  From  these  photo¬ 
graphs,  it  is  possible  to  see  how  shock  wave  systems  are  set  up 
under  various  different  operational  situations. 

As  far  as  turbine  blade  cascades  are  concerned,  in  operating 
conditions  involving  high  subsonic  speeds,  one  will  get  the  produc¬ 
tion  of  localized  areas  of  supersonic  speeds  inside  of  the  passages. 
These  sonic  speeds  first  reach  a  place  somewhere  on  the  back  parts 
of  the  blades  and  then  form  a  boundary  line  on  which  M  =  1  within 
the  passages.  The  lower  reaches  of  this  sort  of  line  reach  from 
the  convex  surface  of  the  backs  of  the  blades  and  product  a  series 
of  expansion  waves  which  cause  the  speed  to  be  increased  somewhat 
more.  This  creates  an  area  in  which  M  >  1.  Within  this  area,  the 
extraneous  reflections  which  are  cast  unto  the  boundary  plane  M  =  1 
by  the  expansion  waves  return  to  the  solid  wall  of  the  back  of  the 
blades  as  compression  waves.  These  compression  waves,  after  they 
make  homogeneous  reflections  on  the  solid  wall  surface,  form  type  \ 
shock  waves.  Because  of  this  fact,  this  localized  area  of  supersonic 
speeds  takes  M  =  1  as  its  front  boundary  and  the  other  boundary  of 
this  small  area — in  the  rear — is  the  \  type  shock  wave  itself. 

This  type  of  phenomenon  is  represented  quite  clearly  and  precisely 
even  in  Schlieren  photographs  in  which  M  is  approximately  0,9. 

Figure  2  is  a  model  of  this  type  of  flow. 


Figure  2.  A  flow  model  of 
localized  areas  of  supersonic 
speeds  under  subsonic  opera¬ 
tional  conditions 

Key:  1 — expansion  wave; 

2 —  compression  wave; 

3 —  boundary  line 


Figure  3-  A  diagram  of  the  cham¬ 
fer  section  of  a  wave  system  in 
a  transonic  planar  blade  cas¬ 
cade  under  supersonic  operating 
conditions . 

Key:  1 — expansion  wave  system; 

2 —  compression  wave  system; 

3 —  wake 


After  the  gas  flow  passes  through  such  localized  areas  of 
supersonic  flow,  it  continues  to  move  along  the  passageways  and 
in  the  vicinity  of  the  exit,  due  to  the  curvature  of  the  trailing 
edges  of  the  air  flow,  it  increases  speed.  In  such  a  situation, 
it  is  easy  (particularly  on  the  backs  of  the  blades  in  the  exhaust) 
to  produce  localized  areas  of  supersonic  speeds. 

When  the  pressure  ratio  behind  the  pressure  in  front  of  the 
blade  cascade  being  considered  and  the  pressure  behind  it  drops  to 
the  critical  level  for  that  ratio,  then  the  chamber-  portion  of  the 
passageways  involved  will  be  formed  into  a  supersonic  wave  system. 
The  model  for  the  flow  In  this  system  is  shown  in  Figure  2. 

The  wave  system  Inside  the  passageways  and  the  exhaust  flow 
field  is  made  up  of  the  original  expansion  wave,  E^ ,  the  reflected 
expansion  wave,  E^,  the  Interior  wake  wave,  K^,  the  reflected  wave 
of  the  interior  wake  wave,  ,  the  exterior  wake  wave,  K^,  as  well 


of  the  blades  Involved. 


as  the  wakes  from  the  trailing;  edges 
After  the  back  pressure  drops  below  the  critical  pressure  for  that 
quantity,  the  throat  section  of  the  passageway  involved  will  reach 
the  speed  of  sound,  the  pressure  at  the  end  point  of  the  interior 
arc.  A,  will  be  the  critical  value  and  in  order  to  reach  a  certain 
back  pressure,  it  is  necessary  to  produce  at  that  point  a  set  of 
expansion  waves,  E^.  This  type  of  expansion  wave  extends  directly 
to  the  back  arc  of  the  neighboring  blades  and  from  there  is  reflected 
out  as  the  set  of  expansion  waves,  E^-  The  gas  flow,  from  the 
back  arc  and  interior  arc  of  the  blade  surfaces  involved,  bends  out 
around  the  trailing  edges.  Moreover,  at  a  certain  point  B  at  a 
given  distance  from  the  trailing  edges,  these  flows  will  come 
together.  Due  to  the  fact  that  the  gas  flow  takes  a  violent  turn 
in  the  vicinity  of  this  point  B,  it  simultaneously  produces  two 
sets  of  compression  wave  systems,  composed  of  a  swallow  tail  wave 
system,  that  is  to  say,  the  interior  wake  wave,  K.^  and  the  exterior 
wake  wave,  K^.  The  interior  wake  wave  has  an  effect  on  the  point  C 
on  the  back  of  the  blade  form  and  mutually  interacts  with  the 
boundary  layers  involved.  After  this  interaction,  the  shock  wave, 

K^,  is  reflected  back  out.  This  shock  wave  is  a  function  of  the 
increase  in  the  M  number  at  the  exhaust.  Concerning  the  increase 
in  strength  of  and  E^ ,  K^,  at  the  point  where  the  jet  enters  the 

back  of  the  blade  form,  moves  in  the  direction  of  the  exhaust  pro¬ 
file.  When  the  expansion  wave  system  puts  its  last  wave  in  a 
basically  parallel  condition  with  the  cross-section  of  the  exhaust, 
then  it  is  possible  to  recognize  the  fact  that  we  have  reached  the 
extreme  limit  of  the  blade  cascade  involved. 

The  reflection  which  takes  place  from  a  wave  system  which  is 
reflected  off  a  solid  wall  is  a  reflection  which  retains  the  form 
of  the  original,  that  is  to  say,  that  the  reflection  of  an  expan¬ 
sion  wave  is  an  expansion  wave.  After  the  reflection  of  a  shock 
wave,  the  reflected  wave  will  still  also  be  a  shock  wave.  This 
point  is  not  only  capable  of  being  verified  from  the  principles  of 
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aerodynamics  [5,6];  it  is  also  completely  verified  by  the  Schlieren 
photographs  which  we  took.  In  the  flow  field  of  the  exhaust,  two 
sets  of  shock  waves  usually  intersect.  If  the  two  sets  which 
intersect  are  compression  waves  of  a  similar  type,  then  after 
they  intersect,  they  come  together  to  form  one  shock  wave  set. 

This  type  of  situation  is  shown  in  its  most  standard  form  in  Plate 
1,  where  the  reflected  wave,  ,  which  is  the  result  of  the  interior 
wake  wave,  K^,  being  in  contact  with  the  solid  wall  surface,  and  the 
exterior  wake  wave,  K^,  come  together  (see  Figure  4a).  The  situa¬ 
tion  in  which  shock  waves  of  different  types  intersect  is  differ¬ 
ent.  Generally  speaking,  the  common  situation  is  one  In  which  the 
two  shock  waves  involved  just  penetrate  each  other.  If  the  strength 
of  the  two  shock  waves  is  equal,  then,  after  the  penetration,  they 
form  a  simple  crossed-fork  pattern,  which  is  a  symmetrical  pattern. 

If  the  strengths  of  the  two  are  not  equal  or  their  M  number  is 
relatively  small,  then  the  vicinity  in  which  these  two  sets  of 
shock  waves  intersect  will  be  turned  into  a  small  section  of  straight 
shock  wave.  The  whole  wave  system  will  present  the  appearance  of 
a  bridge-shape.  It  is  worthwhile  to  point  out  that,  due  to  the 
fact  that  shock  waves  of  different  strengths  are  being  crossed,  the 
gains  in  enthalpy  are  not  the  same  either.  When  the  gas  flow  pene¬ 
trates  two  sets  of  waves,  the  speeds  of  the  top  and  bottom  two 
halves  will  be  different  (although,  at  the  same  time,  the  pressures 
and  directions  of  flow  for  the  two  parts  will  still  be  the  samel. 
Because  of  these  facts,  there  will  exist  in  the  center  a  non-conti- 
nuous  vortical  surface  (see  Figure  4b).  The  gas  flow  model  for  an 
intersection  between  shock  waves  is  as  shown  in  Figure  4.  The  one 
on  the  right  side  is  the  Schlieren  photograph  which  corresponds  to 
this.  On  the  basis  of  the  discussion  above  and  the  analysis  pre¬ 
sented  there,  it  is  possible  to  easily  understand  Plate  1,  which  is 
a  set  of  Schlieren  photographs  which  present  wave  systems  under 
different  operational  conditions. 

2 .  The  form  of  reflected  shock  waves  and  the  Influence  of 
their  convection  flow  characteristics 
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Figure  4.  An  aerodynamic  model  for  the  intersection  of  shock  waves 
and  the  Schlieren  phtographs  for  it. 

Key:  1 — the  intersection  and  remerging  of  sloped  shock  waves  of 

a  similar  type;  2 — the  intersection  of  sloped  shock  waves 
of  different  types;  3 — when  the  strength  of  shock  waves  is 
not  the  same  or  the  M  number  is  relatively  small;  4 — when 
the  strengths  of  the  shock  waves  are  equal. 

Concerning  blade  cascade  wind  tunnels,  due  to  the  fact  that 
the  number  of  blades  is  very  limited,  the  exhaust  flow  field  must 
contain  two  free  boundary  lines.  It  can  be  clearly  shown  from  the 
Schlieren  photographs  which  we  took  that  in  the  reflection  of  a 
shock  wave  from  a  free  boundary,  there  is  a  new  shock  wave  produced. 
This  phenomenon  was  noticed  early  on  by  NACA  [73  and  research  was 
done  on  it  in  1968  by  Von  Karman.  C.  Sieverding  [2]  has  also  given 
his  attention  to  this  problem. 


Concerning  the  phenomenon  discussed  above,  the  experimental 
diagram  which  we  produced  can  be  explained  as  follows.  When  the 
exterior  wake  shock  wave  gets  close  to  the  free  boundary  surface, 
one  sees  formed  an  area  of  transition  from  supersonic  to  subsonic 
speed.  Inside  of  this  area,  the  M  number  gets  gradually  smaller; 
moreover,  this  slowing  is  accompanied  by  the  production  of  vortical 
surfaces  and  an  increase  in  enthalpy.  As  a  result  of  this,  one  sees 
the  shock  wave  involved  gradually  exhibit  curvature.  In  order  to 
maintain  the  equilibrium  of  the  flow,  within  the  transition  area 
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close  to  the  exterior  wake  chock  wave  curvature,  there  will  he 
produced  a  new  shock  wave.  The  Mach  number  reflection  limits 
for  this  model  and  the  classic  model  are  both  similar.  The  latter 
usually  has  a  Mach  number  of  flow  relatively  close  to  1  and  occurs 
when  the  angle  of  folding  is  relatively  large.  Of  course,  in 
order  to  maintain  the  condition  of  equal  pressure  on  the  free 
boundary  surface,  in  the  lower  part  of  the  entering  shock  wave, 
there  is  also  a  set  of  expansion  shock  waves  sent  out.  It  can  be 
seen  from  the  Schlieren  photography  in  plate  1  that: 

1.  In  the  range  of  M}(  —  1  ~  1.1  t  the  wake  area  produced  by 
the  exterior  wake  shock  wave,  K^,  going  through  the  blade  wake  area 
gives  rise  to  a  weak  reflected  shock  wave.  In  response  to  a  reduc¬ 
tion  in  the  ratio  of  the  pressure  in  front  of  a  blade  cascade  and 
the  pressure  behind  it,  the  speeds  in  evidence  at  various  points 

in  the  flow  field  behind  the  exhaust  of  the  blade  cascade  will  be 
increased;  then  the  point  at  which  the  reflected  shock  wave  is  pro¬ 
duced  will  gradually  get  closer  to  the  free  boundary  surface. 

2.  These  reflected  shock  waves  always  return  to  sections  of 
blade  cascade  passage  but  some  sections  of  these  passages  still  do 
not  have  this  type  of  wave.  Because  of  this,  the  periodicity  of 
the  passageways  is  disrupted. 

3.  Following  an  increase  in  the  M  number  at  the  exhaust,  the 
direction  of  the  exterior  wake  shock  wave  manifests  a  change.  The 
direction  of  its  reflected  shock  wave  also  changes  in  response  to 
this  and  the  more  the  process  continues  the  more  the  direction  of 
the  reflected  shock  wave  becomes  parallel  to  the  profile  line  of 
the  exhaust.  The  passageways  which  feel  the  influence  of  the  re¬ 
flected  shock  waves  also  decrease  in  size  in  response  to  these 
other  changes. 

Concerning  the  existence  of  free  surface  reflected  shock  waves, 
these  shock  waves  interfere  with  the  passageways  of  the  blade 
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Figure  5.  Experimental  Curves  for  Transonic  Turbine 
Blade  Cascade  Performance. 
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cascades  involved  and  with  the  periodicity  of  the  exhaust  flow 
field  and  this  creates  a  false  appearance  in  some  areas  which 
har  sere  influence  on  the  reliability  of  the  blade  cascade  wind 
tunnel  data  from  time  to  time.  Figure  5  is  a  display  of  the  aero¬ 
dynamic  characteristics  of  a  certain  blade  cascade.  In  it,  a)  is 
the  relationship  according  to  which  the  loss  coefficient  £  changes 
as  a  function  of  f'^t'  same  diagram,  b)  is  a  distribution 

diagram  of  total  pressure  as  actually  measured  at  various  distances 
along  the  blade  cascade  for  increments  of  M2t *  In  can  be  seen  ^rom 
b)  in  this  diagram  that  the  cascade  losses  are  composed  of  two 
parts.  From  left  to  right,  the  first  loss  area  is  that  formed  by 

« 

the  wake.  The  second  loss  increase  area  is  the  result  of  the 
influence  of  the  wake  edge  shock  wave.  If  we  consult  plate  l's 
Schlieren  photographs,  when  *  1.3,  there  is  a  mutual  influence 
between  the  reflected  shock  wave  from  the  free  boundary  and  the 
exterior  wake  shock  wave  of  blade  cascade  II.  This  effect  causes 
an  enlargement  of  the  second  loss  area.  Because  of  this  fact,  there 
is  an  obvious  increase  in  losses  (see  Figure  5a).  When  M2t  is 

j  •  Is  raised  to  1-35,  the  reflected  shock  wave  moves  toward  the  outside 

of  the  blade  cascade,  leaving  the  passageways,  and  causing  losses 
instead  of  drops.  When  M2t  is  raised  to  reach  1.4,  the  reflected 
shock  wave  reaches  the  back  portion  of  blade  II  which  is  being 
measured  causing  localized  boundary  layer  separation.  There  is  a 
clear  enlargement  of  the  wake  loss  area  and  the  amount  of  loss 
is  again  increased.  When  the  process  of  increase  is  continued 
straight  on  until  M2t  =  1.47  or  thereabouts,  the  reflected  shock 
wave  Involved  moves  out  of  the  back 

portion  of  the  blade  being  measured  and  the  losses  drop  again. 

After  this,  .following  an  increase  in  the  M2t  number,  the  reflected 
shock  wave  reaches  a  direction  which  tends  to  be  parallel  to  the 
profile  line  of  the  blade  cascade  being  considered.  The  influence 
on  the  performance  of  the  passageways  tends  to  disappear.  However, 
due  to  an  increase  in  the  strength  of  the  shock  waves,  the  interior 
,  wake  shock  wave  again  exerts  its  influence  on  the  subsonic  area  at 

the  trailing  edge  and  there  Is  a  sharp  increase  in  the  amount  of 
losses  experienced.  From  the  analysis  advanced  above,  it  can  be 
seen  that  the  loss  curve  for  waves  within  the  range  M.,^  *  1.3  to 
1.47  can  very  possibly  be  a  false  image  produced  by  the  shock  wave 
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reflected  from  the  free  surface .  It  was  also  discovered  during 
these  experiments  that  the  exhaust  angle  of  the  blade  cascades 
being  considered  in  this  section  of  the  operational  ranges  also 
exhibits  wave  movement  and  can  almost  be  explained  by  the  use  of 
this  reason. 

In  order  to  raise  the  reliability  of  the  data,  speaking  in 
terms  of  basics,  one  needs  to  find  a  method  for  eliminating  this 
reflected  shock  wave  phenomenon.  Quite  a  few  scholars  have  put  a 
lot  of  effort  into  this  problem,  coming  up  with  such  solutions  as 
making  use  of  perforated  walls  with  or  without  suction,  using 
movable  final  blades  and  other  such  measures.  However,  the  results 
from  these  attempts  have  not  been  too  impressive.  Looking  at  the 
problem  from  another  angle,  if  one  were  able  to  create  a  mechanism 
by  which  to  carry  out  in-depth  studies  of  reflected  shock  waves, 
then  from  this  basis  it  might  be  possible  to  carry  out  calculations 
of  and  adjustments  for  the  influences  of  blade  cascade  aerodynamic 
characteristics  as  modified  by  reflected  shock  waves.  This  would 
also  make  it  possible  to  some  extent  to  eliminate  the  deficiencies 
of  transonic  planar  blade  cascade  wind  tunnels.  In  fact,  because 
of  the  fact  that  the  position  which  corresponds  to  the  measurement 
point  on  the  shock  wave  involved  is  different  for  each  blade  cas¬ 
cade,  the  influence  which  the  shock  waves  will  have  on  the  con¬ 
struction  of  different  blade  cascades  will  also  certainly  not  be 
limited  to  only  one  kind.  In  our  opinion,  before  we  will  be  able 
to  clearly  understand  the  laws  which  govern  this  type  of  phenomena 
and  before  we  make  precise  experimental  determinations  of  the  char¬ 
acteristics  of  each  different  type  of  blade  cascade  as  far  as  this 
problem  Is  concerned,  it  would  be  very  useful  to  the  objective  of 
raising  the  reliability  of  our  data  to  continue  on  resolutely  with 
the  taking  of  Schlieren  photographs  in  order  to  see  clearly  the 
locations  and  functioning  of  reflected  shock  waves  and  exert  every 
energy  toward  the  choosing  of  appropriate  passageways  and  positions 
for  the  construction  of  blade  cascades. 


3 .  Interference  between  shock  waves  and  boundary  layers 

From  the  Schlieren  photograph  in  plate  1,  it  is  possible  to 
see  that  there  are  several  types  of  interference  phenomena  occur- 
ing  between  shock  waves  and  boundary  layers.  In  the  interval 
when  mu  —  j.o  —  i.2  the  internal  wake  shock  wave  of  the  III  blade 
shoots  into  the  back  arc  of  the  II  blade  and  sends  out  two  shock 
waves.  This  phenomenon  is  the  result  of  interference  between  shock 
waves  and  boundary  layers.  From  (8)  one  can  know  that  the  design 
of  this  set  of  blade  cascades  in  its  constriction  section  is  rela¬ 
tively  good  and  the  gas  flow  along  the  back  arc  of  the  blade  forms 
accelerates  uniformly.  Due  to  the  fact  that  the  acceleration  of  the 
flow  has  a  stabilizing  effect  on  the  boundary  layers  involved,  the 
laminar  flow  flow  configuration  is  capable  of  maintaining  itself 
up  to  a  very  high  critical  Reynolds  number  (Re**  —  l.ox 
Because  of  this  fact,  within  the  range  of  our  experimentation,  the 
boundary  layers  of  the  backs  of  blades  can  still  be  considered  as 
laminar  flows . 

As  is  shown  right  In  Figure  6,  when  the  Interior  wake  shock 
wave  enters  into  the  laminar  flow  of  a  boundary  layer.  It  applies 
to  the  boundary  layer  a  positive  pressure  gradient.  This  type  of 
pressure  gradient,  in  the  subsonic  flow  layer  inside  of  the  bound¬ 
ary  layer,  flows  toward  the  upper  reaches  of  the  flow  to  propagate 
itself.  At  the  places  where  this  shock  wave  enters  the  upper 
reaches,  it  causes  the  beginning  of  an  increase  in  pressure  strength. 
This  type  of  pressure  strength  comes  back  around  and  Is  again  capa¬ 
ble  of  leading  to  an  increase  In  the  thickness  of  the  boundary  layer 
Involved.  Because  of  this  leading  to  localized  changes  in  the  rate 
of  curvature  of  the  flow  lines  involved,  a  system  of  compression 
waves  Is  produced.  This  system  of  compression  waves  follows  behind 
and  comes  together  to  form  the  first  reflected  shock  wave.  It  is 
possible  to  see  that  the  point  of  origin  of  the  first  reflected 
shock  wave  should  be  in  the  upper  reaches  of  the  flow  where  the 
original  point  of  entry  was.  Due  to  the  fact  that  gas  flow  outside 
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Figure  6.  A  model  of  the  interference  which  takes  place 
between  shock  waves  and  boundary  layers  as  well  as  a 
Seblieren  photograph  of  its  typical  form. 

Key:  1 — laminar  flow  boundary  layer;  2 — the  incoming  shock  wave; 

3 — turbulence  flow  boundary  layer;  4 — expansion  wave  system; 
5 — the  first  reflected  shock  wave;  6 — the  second  reflected 
shock  wave 
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Figure  7.  A  flow  model  of  the  interference  between 
shock  waves  and  boundary  layers. 

Key:  1 — normal  reflection;  2 — fork-type  reflection;  3 — turbulence 

flow  boundary  layer;  4 — shock  wave  entry;  5 — reflected 
shock  wave;  6 — area  of  separation 
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of  the  boundary  layer  penetrates  the  shock  wave  at  points  where 
its  pressure  strength  is  larger  than  that  in  the  boundary  layer, 
the  pressure  strength  gradient  which  this  creates  will  again 
cause  a  folding  back  of  the  flow  lines  toward  the  wall  surface 
and  will  lead  to  an  increase  in  the  curvature  of  the  surface  of 
the  blade  being  considered.  From  all  of  this  comes  the  production 
of  a  set  of  expansion  waves  that  constitute  a  system.  However, 
very  frequently,  due  to  the  fact  that  the  angle  of  deflection  of 
the  expansion  waves  is  too  great,  this  tends  to  cause  the  air  flow 
to  move  in  toward  the  wall  surface  and  the  wall  surface  further 
down  has  an  effect  of  adjusting  the  direction  of  the  gas  flow.  It 
is  as  if  one  had  a  concave  angle  there  leading  to  the  production 
of  a  system  of  compression  waves  which  later  come  together  to  form 
the  second  reflected  shock  wave.  Following  hard  on  after  the 
second  shock  wave,  the  boundary  layer  follows  along  and  turns  to 
become  the  turbulence  flow.  Figure  6  is  a  model  of  the  interfer¬ 
ence  between  shock  waves  and  boundary  layers  as  well  as  a  Schlieren 
photograph  of  the  typical  form  of  this  phenomenon. 

When  a  shock  wave  makes  its  effects  felt  in  a  turbulence  flow 
boundary  layer,  the  situation  is  very  much  different.  Due  to  the 
fact  that  the  turbulence  flow  boundary  layer  Is  capable  of  rela¬ 
tively  well  overcoming  the  reverse  pressure  gradient,  the  effects 
produced  by  the  Interference  are  much  smaller.  At  this  time  the 
thickness  of  the  boundary  layer  involved  Increases  only  very 
slightly.  In  the  flow  above  and  below  the  point  of  entry,  one  will 
see  the  production  of  several  beams  with  very  small  ranges  of 
influence.  Because  of  all  this,  the  reflection  of  shock  waves  in 
a  turbulence  flow  boundary  layer  closely  approximates  a  normal 
reflection,  that  is  to  say,  the  angle  of  entry  of  the  shock  wave 
and  the  angle  of  reflection  will  be  approximately  equal.  The  model 
for  this  type  of  flow  is  as  shown  in  Figure  7a.  The  corresponding 
Schlieren  photographs  are  as  shown  in  plate  1  for  the  M2t  *  1.3 
'v  1.35  situation. 


When  K2  continuously  increases,  the  strength  of  the  entering 
shock  wave  increases.  The  effect  which  the  entering  shock  wave 
will  have  on  a  turbulence  flow  will  cause  the  occurrence  of  a 
small  amount  of  separation.  The  reflected  shock  wave  which  occurs 
at  this  time  will  begin  from  the  point  of  entry  and  spread  up  the 
flow  forming  a  fork-type  reflected  shock  w ave  as  shown  in  Figure 
7b.  The  Schlieren  photograph  for  this  type  of  reflection  can  be 
seen  in  plate  1  as  the  situation  in  which  M2t  =  1.4. 

Due  to  the  fact  that  the  interference  between  shock  waves  and 
boundary  layers  is  a  type  of  irregular,  complicated  physical  phen¬ 
omenon,  at  present,  in  the  pictures  which  are  taken  of  turbine 
blade  cascades  in  transonic  applications  still  lack  clarity.  When 
we  went  through  the  Schlieren  photographs  for  this  set  of  blade 
cascades,  it  was  only  possible  to  make  a  few  very  crude  observa¬ 
tions  and  come  to  some  crude  understandings  of  these  phenomena. 

The  analysis  and  editing  involved  in  this  article  were  done 
by  Min  Da-fu. 
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THE  SHOCK  WAVE  STRUCTURE  AND  THE  ANALYSIS  OF 
THE  FLOW  THROUGH  TRANSONIC  TURBINE  CASCADE 

Min  Da-fu 

(The  mir  Lab.  0/  the  Beijing  Heavily  Electrical  Machine  Works) 

v 

Abstract 

Throughout  a  let  of  Schlieren  photos  from  the  aerodynamic  experiments  of  the 
Transonic  Tunbine  cascade,  the  mechanism  of  the  generation  and  development  of  the 
shock  waves  in  channel  and  exit  flow  field  of  the  cascade  is  studied.  In  these  photos, 
there  are  a  series  of  deflect  shock  waves  from  the  free  boundary  of  the  jet,  Passing 
through  the  analysis  of  this  physical  phenomenon,  the  satisfied  explanation  on  the 
aerodynamic  performance  which  hus  a  multipeak  type  of  curve  is  obtained.  In  the  3rd 
of  this  article,  three  types  of  interaction  between  shock  waves  and  boundary  layer  are 
observed  and  discussed  simply. 
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A  THEORETICAL  ANALYSIS  OF  THE  STREAK  SURFACE  OF  REVOLUTION 
WITH  SUPERSONIC  INLET  FLOW  IN  A  TRANSONIC  AXIAL  COMPRESSOR* 

Xu  Jian-zhong,  Jiang  Zuo-ren,  Yang  Jin-sheng,  Zhang  Ying  and 
Du  Zhuan-she 

(Institute  of  Engineering  Thermophysics,  Academia  Sinica) 

If  we  start  with  changes  in  the  direction  and  numerical  value 
of  the  speeds  corresponding  to  gas  flows  when  going  through  a  shock 
wave  as  these  changes  lead  to  jumps  in  the  absolute  speeds  of  gas 
flows  as  well  as  to  jumps  in  the  tangential  components  of  these  gas 
flows,  one  should  make  use  of  planar  oblique  wave  relationships  to 
establish  a  formula  for  expressing  A(Vg  ),  losses  and  total  press¬ 
ure  ratios  and  from  the  applicable  theory,  we  explained  the  impor¬ 
tant  phenomena  related  to  gas  pressure  increases  when  passing 
through  shock  waves.  We  considered  the  results  of  experimentation, 
assumed  the  positions  at  which  shock  waves  were  crossed  as  well  as 
their  configurations  and  made  use  of  the  characteristic  curve 
method  in  the  intake  area.  We  also  made  a  theoretical  analysis  of 
the  surface  of  revolution  of  the  supersonic  intake.  From  this  anal¬ 
ysis  it  could  be  seen  that  there  were  several  special  points  asso¬ 
ciated  with  transonic  compression  and  compressors.  Calculations 
also  showed  that  in  the  center  line  exhaust,  when  the  inverse  pro¬ 
blems  of  the  velocity  moment  Vgr  and  the  S2  flow  are  calculated,  the 
numerical  values  which  are  arrived  at  for  the  two  are  very  close. 

If  one  arranges,  for  the  blades  Involved,  the  average  results  of 
calculating  regular  values  of  S2  and  the  Mach  number  distributions 
which  are  obtained  by  substituting  into  the  S2  flow  surface  calcula¬ 
tions  the  flow  filament  thickness  and  the  *V»  which  are  gotten 
by  considering  the  shock  waves  involved,  then  there  will  be  a  parti¬ 
cularly  large  difference  between  conditions  in  front  of  the  shock 
waves  involved  and  the  conditions  behind  them. 

— 

This  article  was  read  at  the  3rd  Annual  All-China  Engineering 

Thermophysics  and  Technology  Conference  held  in  Gui  Lin  in  April 

of  1980. 
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1 .  FOREWORD 

In  order  to  make  an  in-depth  study  of  the  complex  three- 
dimensional  flows  inside  of  transonic  axial-type  compressors,  we 
took  its  aerodynamic  and  thermodynamic  design  as  offering  a  reli¬ 
able  theoretical  foundation  and  on  the  basis  of  the  theory  and 
methods  presented  in  [1],  we  designed  a  single-stage  research  com¬ 
pressor  [2]. 

As  far  as  transonic  flows  which  carry  shock  waves  with  them 
are  concerned,  there  are  two  methods  of  understanding  them.  One 
type  of  method  is  to  do  numerical  calculations  and  obtain  the  whole 
flow  field  after  one  has  established  a  physical  model  and  set  up 
the  necessary  mathematical  equations  as  well  as  establishing  the 
conditions  for  the  solutions  desired.  When  this  type  of  solution 
is  employed,  the  jumps  in  the  parameters  behind  the  wave  are  re¬ 
placed  by  continuous  changes  and  the  shock  wave  is  pulled  apart 
into  two  or  three  coordinate  networks.  The  other  method  is  to  assume 
the  position  and  configuration  of  the  shock  waves  involved  and  solve 
for  the  flow  fields  in  front  and  behind  them.  The  calculations 
involved  in  this  type  of  solution  are  relatively  simple.  The  key 
to  using  this  method  is  the  accuracy  of  the  determinations  for  the 
position  and  configurations  of  the  shock  waves  involved  and  these 
quantities  need  to  be  determined  by  experimentation. 

References  [3*^]  make  use  of  the  laser  method  for  measuring 
speed.  These  references  carry  out  relatively  detailed  measurements 
of  the  flow  fields  of  transonic  axial  compressors  which  contain 
shock  waves  in  them  and  make  a  comparison  between  these  results  and 
the  results  obtained  by  normal  measurement  and  flow  calculation 
techniques.  These  references  obtained  quite  a  few  interesting 
results . 

In  order  to  analyze  the  effect  which  shock  wave  have  on  gas 
flow,  this  article  starts  with  changes  in  the  gas  flow  velocity 


36 


triangle  behind  the  shock  wave  being  considered,  then  use  is  made 
of  the  planar  oblique  relationship  formula  and  we  deduce  the  cor¬ 
responding  values  of  a(Vq„)  ,  pressure  losses  and  total  nr^ssure  ratios 
which  are  expressed  by  this  type  of  formula.  This  explains  the 
important  phenomenon  involved  when  the  supersonic  intake  cross- 
section  shock  waves  of  transonic  compressors  cause  additional  press¬ 
urization  of  gas  flow  or  "supercharging".  On  this  foundation,  we 
consult  the  results  of  experimentation,  assume  certain  positions 
and  configurations  for  the  shock  waves  involved.  We  then  carry 
out  aerodynamic  and  thermodynamic  analyses  and  calculations  on  the 
top  cross-sections  of  the  rotors  in  [2]  and  then  discuss  the  results 
obtained. 

2.  CHARACTERISTIC  CURVE  CALCULATIONS  FOR  THE  SHOCK  WAVES  IN  FRONT 
OF  BLADE  CASCADES  AND  FOR  THE  INTAKE  AREAS  OF  THESE  CASCADFS 

In  order  to  carry  out  aerodynamic  and  thermodynamic  calcula¬ 
tions  of  the  intake  surface  of  revolution  and  supersonic  condi¬ 
tions,  it  is  necessary  to  accurately  know  the  position  and  config¬ 
uration  of  the  shock  wave  in  front  of  the  blade  cascade  involved. 

In  this  regard,  as  far  as  Isolated  wing  forms  are  concerned,  we 
already  have,  on  the  basis  of  experimental  data,  the  relatively 
easily  adaptable  MOECKEL  model  [5].  Obviously,  as  far  as  blade  cas¬ 
cades  in  compressors  are  concerned,  one  should  give  consideration 
to  the  lack  of  uniformity  in  the  flow  fields  of  the  upper  reaches 
of  the  flows  involved.  Other  such  factors  should  also  be  similarly 
considered  in  the  light  of  their  influence  on  the  make-up  of  the 
position  and  configuration  of  shock  waves  in  front  of  blade  cascades. 
These  other  factors  are  such  things  as  the  mutual  Interaction  between 
blades  and  the  structure  of  the  blade  cascade  itself. 

Concerning  a  rotational  supersonic  intake  area  from  behind  a 
detached  shock  wave  to  the  passageways  in  front  of  the  shock  waves, 
it  is  possible  to  make  use  of  the  characteristic  curve  method  in 
order  to  solve  for  the  quantities  desired.  Reference  [6]  gives  us 
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a  characteristic  curve  equation  which  includes  the  flow  filament 
thickness  t  and  the  redius  r  of  the  surface  of  revolution  as  well 
as  other  relationships  on  which  this  equation  is  based,  that  is 
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Based  on  these  equations,  starting  out  from  an  initial  line 
which  is  outside  of  the  sonic  line  of  the  detached  wave  involved, 
and  on  the  basis  of  the  corresponding  boundary  conditions,  it  is 
possible  to  solve  for  the  supersonic  Intake  flow  field. 

We  know  that  when  the  speed  which  corresponds  to  intake  is 
supersonic  and  the  meridional  component  of  it  is  subsonic,  it  is 
possible  for  the  disturbances  produced  in  the  blades  to  spread  in 
the  direction  of  the  upper  reaches  of  the  flow.  Because  of  this. 

In  front  of  the  blade  cascade  one  sees  set  up  a  condition  of  stable 
periodicity.  In  this  way,  using  the  characteristic  curve  method, 
beginning  from  the  assumption  of  a  semi-infinite  blade  cascade,  we 
went  through  the  calculations  Involved  for  three  or  four  passageways. 
In  these  flow  fields  the  M  numbers  and  the  B  values  for  the  corres¬ 
ponding  positions  were  all  the  same  and  satisfied  the  periodicity 
conditions.  Moreover,  when  we  included  in  the  range  a  place  which 
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was  infinitely  distant  from  the  Liade  cascade,  the  tangential  dis¬ 
tance  was  one  cascade  length  and  this  type  of  control  body  was 
set  up  as  the  basis  for  the  basic  aerodynamic  and  thermodynamic 
equations  necessary  in  order  to  solve  for  even  and  uniform  condi¬ 
tions  of  incoming  flow.  After  carrying  out  quite  a  number  of  this 
type  of  calculations,  it  became  possible  to  precisely  determine 
unique  angles  of  intake.  From  this,  it  is  not  difficult  to  obtain 
blade  forms  which  satisfy  unique  angle  of  attack  relationships  and 
at  the  same  time,  we  also  solved  for  the  flow  field  of  the  super¬ 
sonic  intake  area. 

3.  PASSAGEWAY  SHOCK  WAVES 


Normally,  passageway  shock  waves  are  extensions  of  shock  waves 
in  front  of  the  blade  cascade  in  the  direction  of  blade  cascade 
passageways  inside.  As  far  as  such  conditions  of  these  passageway 
shock  waves  as  position,  configuration  and  strength  are  concerned, 
we  have  carried  out  a  good  deal  of  experimental  research  in  this 
area.  On  the  basis  of  the  results  in  [3,^3,  generally  speaking, 
what  appears  in  the  inside  of  passageways  are  oblique  shock  waves. 

In  order  to  accurately  analyze  the  situation  in  which  we  find 
the  gas  flow  in  front  of  and  behind  the  shock  wave  being  considered, 
strictly  speaking,  because  of  the  qualities  of  space  involved,  one 
should  make  use  of  the  shock  wave  relationships  presented  in  [73. 

In  order  to  obtain  formulae  for  expressing  problems  that  are  easy 
to  explain,  this  article  still  makes  use  of  simple  relationships  of 
the  planar  oblique  type. 

As  far  as  planar  oblique  shock  waves  are  concerned,  there  are 
two  types  of  different  situations.  One  of  these  is  the  .oblique 
shock  wave  in  which  B2  =  B-j_  -  e  is  the  case  (.Figure  1).  Besides 
this,  e  is  the  angle  of  folding  for  the  corresponding  air  flow.  It 
is  determined  from  the  air  flow  parameters  for  the  area  in  front  of 
the  waves  concerned,  and  is 


39 


ctge  — 


A/?  *inJ  f4y 


(l-l’XMidn’/H-l) 


—  'g 


(1) 


In  this  equation  A1  -  U  —  0/(4  +  0.  represents  the  included 

angle  between  the  incoming  flow  and  the  shock  wave  involved. 

After  going  through  the  shock  wave,  the  magnitude  of  the  speed 
of  the  gas  flow  also  exhibits  changes,  that  is 
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From  the  velocity  triangle  which  is  shown  in  Figure  1,  one 
can  see  that  when  going  through  shock  waves,  changes  in  the  numer¬ 
ical  values  for  the  speeds  of  the  corresponding  gases  as  well  as  in 
their  direction  of  flow  cause  changes  in  the  absolute  numerical 
velocity  values  for  these  gases  as  well  as  in  their  directions,  and 
one  can  see  the  production  of  jumps  in  the  tangential  components 
concerned.  From  the  figure  one  can  learn  that 
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.  „  f.  sin  ft  1  +  *2(A/j  sinVi  —  1) 
smft  II  — 


—  sin | 


sin  ft 


M }  sin  f*\  sin  f*i 


] 


—  2rlF,slnft(I  +  ctg  ft  ctg  f«i)(  sin2pi  —  l/A/{)/(^  + 


(2) 


From  these  equations,  one  can  see  that  the  jumps  in  the  velo¬ 
city  moment  of  the  gas  flow  being  considered,  a(i>)  and  A/i, « 
are  related  to  0^.  Most  particularly  in  the  case  of  the  positive 
shock  wave  In  the  illustration,  w  —  »r/2 ,  ft  ft  “  ^  ,  and  on  this 

basis,  we  can  obtain 

A( IV )  -  2rW'liiin0(l  -  1/Af})/U  +0  (3) 


When  passing  through  this  positive  shock  wave,  even  though  the 
corresponding  gas  flow  does  not  give  rise  to  bending,  due  to  jumps 
In  the  numerical  values  for  the  velocities  involved,  it  does  give 
rise  to  jumps  in  the  value  of  £(*V)  . 
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Another  type  of  situation  is  the  one  shown  in  Figure  2,  that 
is,  &  —  /?,  +  £  .  In  such  a  situation,  A(tV)  is  similarly 

related  to  M^,  and  6^,  that  is 

A(  V ur)  —  2rW i  sin  /?j(  1  —  Ctg  ctg  /i,)(  ^in2  /j,  —  l/A/{)/({  +  O  (  4  ) 

What  we  are  primarily  concerned  with  in  this  type  of  situation 
as  far  as  the  orthogonal  shock  wave  is  concerned,  is  the  simplifi¬ 
cation  of  (4)  into  (3). 


As  far  as  these  two  types  of  oblique  shock  waves  are  concerned, 
the  absolute  total  pressure  ratios  of  the  gas  flow  in  front  of  the 
wave  and  the  gas  flow  behind  the  wave  can  both  be  represented  as 


Pi* 

Pic 


*±L/ 

(*+  !)*-»( 


_ 1 _ 

2*Ai}sin1/n-*  +  1 


* 


X  ■ 

Mj  sin2  ft, 

{  —  1  to  I 

l+  *  RTlo  ’ 

(-1 

;  —  1 )  A/  {  sin2  ft,  +  2 

(5) 


In  this  equation,  the  quantity  A(*V)  can  be  calculated  by 
dividing  up  the  task  between  Equations  (2)  and  (4). 

It  can  be  seen  that  the  larger  the  M  number  in  front  of  the 
wave  is,  the  larger  a(JV)  Is,  and  the  larger  the  ratio  of  total 
absolute  pressures  P2c/P10  is‘ 


If  we  make  use  of  the  corresponding  total  pressure  ratio  of 
gas  flows  in  order  to  represent  losses  produced  after  going  through 
an  oblique  shock  wave,  then  It  is  possible  to  write 


Ko 

P',0 


1)* 


/ _ 1 _ M\  sin*  m 
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Figure  1.  The  oblique  shock 
wave  32=Bi-e  and  the  velocity 
triangle  for  the  gas  flow  both 
in  front  of  and  behind  the 
wave 

Key:  1 — shock  wave 


Figure  2.  The  oblique  shock 
wave  3z=6i+e  and  the  velocity 
triangle  for  the  gas  flow  both 
in  front  of  and  behind  the 
wave 


One  obvious  characteristic  of  this  equation  is  that  the  losses 
are  only  determined  by  sin  y-^,  and  are  not  related  to  In 

this  way,  we  are  easily  capable  of  going  through  aerodynamic  designs 
in  order  to  produce  different  kinds  of  shock  wave  systems  so  that 
the  losses  produced  are  relatively  small  and  A(_VQ  )  is  relatively 

O  I* 

large.  Of  course,  the  losses  we  are  talking  about  here  are  still 
losses  which  are  produced  by  going  through  the  body  of  a  shock  wave 
itself.  However,  it  can  be  hoped  that  the  separation  losses  which 
are  caused  by  shock  waves  which  are  relatively  weak  will  also  be 
relatively  small. 
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Figure  3.  Flow  Path  Diagram  Figure  4,  Distribution  of  Svv  Along 

the  Passage  Centerline  on  Surface 
of  Revolution 
1  -  surface 

From  (5)  and  (.6),  it  is  possible  to  obtain  the  relationship 
set  out  below  which  connects  the  absolute  total  pressure  ratios 
for  gases  in  front  of  and  behind  a  shock  wave  to  the  correspond¬ 
ing  total  pressure  ratios,  that  is 


Pio 


What  is  interesting  is  the  fact  that  even  though  Pio/Ko  < 1  , 
Pjo/Pio>1  is  still  the  case. 

4.  CALCULATION  RESULTS  AND  DISCUSSIONS 

Concerning  the  use  of  the  characteristic  curve  method  in  the 
Intake  area,  one  postulates  a  position  and  configuration  for  a 
passageway  shock  wave  and  then  makes  use  of  a  unique  flow  assump¬ 
tion  which  takes  into  consideration  the  flow  thickness  and  losses 
Involved  and  applies  this  to  the  gas  flow  behind  the  wave.  Con¬ 
cerning  the  carrying  out  of  calculations  of  the  II  cross-sections 
and  IV  cross-sections  of  the  tips  of  the  rotors  from  [2],  that  is 
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Figure  5-  Changes  in  the  II  sur-  Figure  6.  Changes  in  the  IV 
face  of  revolution  gas  flow  angle  surface  of  revolution  gas  flow 
along  the  center  line  of  a  pass-  angle  along  the  center  line  of 
ageway.  a  passageway. 

Key:  1 — center  arc  line; 

2 — calculated  result 

Figure  3,  we  will  take  the  results  from  these  calculations  and  put 
them  together  with  the  calculation  results  derived  for  the  subsonic 
surface  of  revolution  and  then  take  them  both  and  move  them  into 
the  center  of  the  flow  surface  S2,  that  is  ( s*)  .  Then,  we  will 

carry  out  a  new  set  of  calculations. 

Figure  4  is  a  representation  of  changes  in  the  surface  of 

revolution  II  and  IV  along  a  passage  centerline  V0p.  Figures  5  and 

o  are  respectively  the  changes  In  the  gas  flow  angles  along  the 

center  lines  of  passageways  which  correspond  these  two  surfaces  of 

revolution.  The  Important  point  is  this.  The  VQ  value  for  air 

wr 

flows  which  pass  through  shock  waves  steeply  and  violently  increases. 
For  the  intake  cross-section  IV  with  a  relatively  low  M  number, 

ACV0r)  reaches  99 55  of  the  A(V0r,)  value  for  the  entire  blade  cascade 
passageway  Involved.  In  the  case  of  cross-section  II,  in  which 
VQr  exceeds  the  numerical  value  in  the  exhaust,  causing  the  air 
flow  along  the  center  line  of  the  passageway  behind  the  wave  to 
accelerate,  we  see  the  occurrence  of  a  turbine  situation  in  which 
^0r  again  gradually  decreases  until  it  reaches  the  exhaust  value. 
Correspondingly,  when  passing  through  an  oblique  shock  wave,  the 
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comparison  of  the  results  from  various  different 
methods  used  for  calculating  the  II  cross-section 
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Key:  1 — quantity;  2 — calculation  method;  3 — two-dimensional 

analysis;  4 — orthogonal  shock  wave;  5 — oblique  shock  wave; 

6 — characteristic  curve  oblique  shock  wave;  7 — front;  8 — 
behind;  9 — front;  10 — behind;  11 — wave;  12 — exhaust; 

13 — total 

corresponding  angle  of  air  flow  along  the  center  line  of  the  pass¬ 
ageway  involved  drops  abruptly.  In  the  IV  cross-section,  these 
values  are  somewhat  higher  than  the  numerical  values  in  the  exhaust; 
moreover,  in  the  II  cross-section,  they  are  somewhat  lower  than 
the  ^numerical  values  in  the  exhaust  section;  in  the  passageways 
that  follow,  these  values  gradually  come  back  up  and  reach  the 
values  which  are  found  in  the  exhaust.  Obviously,  if  the  shock 
wave  concerned  is  an  orthogonal  one,  then  the  gas  flow  angle  behind 
the  wave  will  not  change. 

What  is  worthy  of  note  is  the  fact  that  after  we  make  use  of 
the  oblique  shock  wave  model  in  each  of  the  various  cross-sections 
in  the  exhaust  area,  when  power  is  added  up  and  the  calculation  of 
the  inverse  problem  in  S£  is  calculated,  the  numerical  values  which 
one  gets  are  relatively  close  to  each  other,  the  error  being  inside 
a  range  of  3?.  When  we  make  use  of  two-dimensional  analysis  of 
shock  waves  in  our  calculations,  this  error  is  very  large  (Tables 
1  and  2).  This  shows  the  fact  that  the  oblique  shock  wave  model  is 
capable  of  giving  excellent  results  and  the  reason  for  this  is  that 
it  molds  itself  relatively  more  closely  to  the  actual  situation  of 
the  shock  waves  in  the  vicinity  of  the  center  lines  in  the  passageways 
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TABLE  2.  A  comparison  of  the  calculated  results  obtained  by 
the  use  ol’  different  methods  for  cross-section  IV. 
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Figure  7.  The  distribution 
of  the  surface  of  revolution 
1  P  -  t  along  the  center 
5  P  line  of  the  passage¬ 
way  involved. 

Key:  1 — surface  of  revolution  II 
2 — surface  of  revolution  IV 
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of  the  blade  cascades  being  studied.  The  main  reason  is  also  as 
follows:  On  the  center  line  of  the  passageways  involved,  the  gas 
flow  behind  the  shock  waves  concerned  Is  in  a  turbine  configuration 
however,  this  does  not  mean  that  the  entire  passageway  Is  necess¬ 
arily  in  the  same  condition.  In  fact,  in  these  passageways,  the 
various  configurations  and  strengths  of  shock  waves  are  extremely 
different  from  each  other.  Because  of  this  fact,  generally  speak¬ 
ing,  in  order  to  accurately  solve  this  type  of  flow,  one  ought  to 
do  the  sort  of  thing  pointed  out  in  [73  that  is  to  say,  make  simul¬ 
taneous  calculations  of  several  S-^  flow  surfaces  and  several  S£ 
flow  surfaces. 


From  Figure  M ,  it  can  also  be  seen  that,  even  though  the 
intake  M  number  of  the  II  cross-section  is  relatively  high,  it  is 
only  necessary  to  design  the  structure  in  an  appropriate  way  and 
it  becomes  possible  for  the  M  number  in  front  of  the  oblique  shock 
wave  concerned  to  not  be  excessively  high  for  the  losses  of  this 
type  to  remain  small  and  for  one  to  obtain  the  required  total 
pressure  ratios. 

Figure  7  is  a  display  of  the  changes  in  (P  —  /)/£F  for  cross- 
sections  II  and  IV  along  the  center  line.  This  Figure  points  out 
that  this  parameter,  when  it  passes  through  an  oblique  shock  wave 
also  produces  abrupt  changes  in  its  values.  This  interesting  result 
clearly  points  out  the  fact  that  the  flow  filament  S2  behind  the 
shock  wave  concerned  gives  rise  to  deflections. 
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Figure  8.  A  comparison  of  the 
calculated  standard  through  cur¬ 
rent  for  the  surface  of  revolu¬ 
tion  II  and  S2  calculated  Mach 

number.  c 


Normal  Calculation  of  the  S2  Flow 

Figure  9.  A  comparison  of  the 
standard  calculated  through  flow 
for  the  surface  of  revolution  IV 
and  the  calculated  S2c  Mach  number. 


Finally,  Figures  8  and  9  present  the  corresponding  distribu¬ 
tions  of  M  numbers  which  are  obtained  for  cross-sections  II  and  IV 
by  making  use  of  calculations  of  the  standard  through  flow  and  cal¬ 
culations  of  the  S2C  which  takes  into  consideration  oblique  shock 
waves.  It  is  possible  to  see  that  in  the  interval  between  the 
front  and  the  back  of  the  outside  of  blades,  the  errors  in  these 
two  quantities  are  not  large.  However,  inside  the  blades  the  two 
quantities  are  quite  distinct  from  each  other.  This  is  particularly 
true  in  the  area  which  is  in  the  vicinity  of  the  shock  wave.  In 
front  of  the  wave,  the  M  number  which  is  calculated  on  the  basis  of 
the  standard  method  is  much  smaller  when  compared  to  its  companion 
quantity. 

Behind  the  wave,  it  will  be  larger  than  its  companion  quantity. 
It  can  be  seen  from  this  that  in  order  to  cause  the  air  flow  cal¬ 
culations  for  a  transonic  compressor  as  well  as  its  design  closely 
match  the  actual  situation,  it  is  necessary  to  make  use  of  the  dis¬ 
tributions  of  the  quantities  v»r>S  and' (P  —  «)/£Pas  they  were  calculated 
with  due  consideration  given  to  the  shock  waves  involved,  and  use 
these  distributions  to  transfer  into  the  calculation  of  S_  .  More- 

c.  C 

over,  looking  from  the  point  of  view  of  our  calculations,  in  the 
area  in  the  vicinity  of  the  shock  waves  involved,  the  S2C  calcula¬ 
tion  stations  ought  to  be  selected  in  larger  numbers.  In  this 
way,  the  calculations  for  the  area  in  the  vicinity  of  the  shock 
wave  would  become  even  more  accurate. 

5.  CONCLUSIONS 

Passing  through  a  shock  wave,  the  numerical  values  of  velo¬ 
cities  which  correspond  to  gases  under  consideration  as  well  as  the 
direction  of  movement  of  these  gases  exhibit  changes.  These  changes 
lead  to  changes  in  the  absolute  numerical  values  for  velocity  and 
direction  of  movement  as  these  relate  to  the  same  gases.  Because 
of  these  factors,  one  sees  the  occurrence  of  jumps  in  the  tangential 
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velocities,  and  these  jumps  produce  strong  and  violent  compression 
effects  on  the  gases  involved. 

The  results  of  calculations  show  that  on  the  center  lines  of 
passageways  there  are  abrupt  drops  in  the  gas  flow  angles  which 
correspond  to  the  passage  through  the  shock  wave  or  waves  involved. 
These  angles  of  flow  are  sometimes  lower  than  the  angles  of  flow 
at  the  exhaust.  Moreover,  A(Vg^)  exhibits  sudden  increases,  even 
to  the  point  where  it  is  possible  for  its  values  to  exceed  the 
values  found  at  the  exhaust.  The  numerical  values  of  Vgr  at  the 
exhaust  and  the  values  which  come  out  of  calculations  of  S2c  vary 
from  each  other  in  a  range  of  about  3% •  At  the  same  time,  the 
results  of  calculations  of  (P-t)/£P  point  out  that  when  passing 
through  a  shock  wave,  the  flow  surfaces  exhibit  deflection  and 
twisting.  This  is  an  important  point. 

Concerning  the  taking  of  the  results  of  calculating  the  surface 
of  revolution  by  using  this  type  of  oblique  shock  wave  model  and 
substituting  them  into  the  calculations  of  S2c,  the  distribution 
of  parameters  which  one  obtains  in  this  way  exhibit  very  large 
differences  from  the  amounts  of  flow  calculated  as  standard  condi¬ 
tions  and  those  present  on  the  inside  of  the  blades.  This  dis¬ 
crepancy  is  particularly  pronounced  for  values  calculated  in  the 
areas  in  front  of  shock  waves  and  values  calculated  for  areas 
behind  them. 
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AN  THEORETICAL  ANALYSIS  ON  THE  STREAM  SURFACE  OF 
REVOLUTION  WITH  SUPERSONIC  INLET  FLOW  IN  A 
TRANSONIC  AXIAL  COMPRESSOR 


Xu  Jian-zhoug  Jiang  Zuo-ren  Yang  Jin-sltetig  Zlmng  Ying 

Du  Zhutpi-she 

(Institute  of  Engineering  Thermo  phy  ties,  Academia  Sinica) 


Abstract 

Starting  from  jumps  of  absolute  velocity  anti  its  tungential  component  due  to 
changes  of  magnitude  and  direction  of  the  relative  velocity  while  passing  through  a 
shock  aud  using  the  relations  for  a  plane  oblique  shock,  the  expressions  for  A(Vjr) 
loss  and  total  pressure  ratio  of  gas  are  derived  aud  the  fact  that  the  gas  is  supereliarged 
across  the  shock  is  explained  theoretically.  An  uerotliermodynamic  analysis  on  the 
stream  surface  of  revolution  with  supersonic  inlet  flow  is  made  by  assuming  the  location 
and  shape  of  the  pussuge  shock  and  applying  the  characteristics  method  in  the  inlet 
region.  Some  of  important  features  in  a  transonic  compressor  can  be  seen  from  the 
analysis.  The  moment  of  velocity  from  such  calculations  is  close  to  that  given  in  the 
calculation  of  inverse  problem  ulong  a  central  8,  stream  surface.  There  are  great  dif¬ 
ferences  in  Mach  number  distributions  between  the  calculations  of  the  central  S,  stream 
surface  with  the  variations  of  V,r-As/li  and  the  stream  filament  thickness  in  which  the 
shock  relations  arc  taken  into  account  as  its  input,  and  of  the  conveutioual  average  8% 
stream  surface.  1  1 
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PRELIMINARY  RESEARCH  ON  ABLATION  EXPERIMENT  TECHNIQUES  IN 
TURBULENT  FLOW  OVER  A  FLAT  PLATE  SPECIMEN* 

Han  Yin-da,  Wang  Ke-xiang,  Yu  Guang-rong 

Chinese  Aerodynamic  Research  and  Development  Center  (CARDC)) 

SUMMARY 

This  article  involves  the  test  production  of  a  new  type  of 
method  involving  an  electric  arc  heated  supersonic  flow  placed 
at  a  small  angle  to  a  flat  plate  specimen.  This  type  of  method 
is  capable  of  overcoming  the  inadequacies  of  the  pipe  test  method 
which  is  employed  now.  The  speed  of  ablation  which  is  measured 
in  this  method  is  relatively  uniform.  Moreover,  it  is  capable  of 
easily  measuring  the  surface  temperature  and  surface  radiant  heat 
flow.  Also  in  this  article  there  is  a  discussion  of  experimental 
principles,  the  structure  of  experimental  equipment,  measurement 
parameters  and  calculation  methods.  This  study  was  completed  in  a 
measurement  arc  chamber  in  which  pressures  ranged  from  10.2  to 
17.3  kg/cm2  and  enthalpies  of  from  520  to  2 h 0  kcal/kg.  This  arti¬ 
cle  gives  the  results  of  ablation  experiments  with  teflon  models 
under  five  types  of  turbulent  conditions. 


This  article  was  read  at  the  2nd  Annual  All-China  Engineering 
Thermal  Physics  Conference  in  Hang  Zhou  in  November  of  1978. 
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SYMBOLS 


a 

sonic  speed;  entbalpic  loss 
coefficient  behind  the 

P  : 

density 

throat  section  of  a  jet 

u  : 

viscosity  coefficient 

tube 

a  : 

air  constituent 

Ga 

air  flow  amount 

t  : 

shear  stress 

Gw 

amount  of  cooling  water  flow  e: 

full  radiation  coefficient 

H 

0 

total  enthalpy 

a  : 

Boltzmann  constant 

Hr 

recovery  enthalpy 

Superscripts : 

h 

static  enthalpy 

«  . 

reference  enthalpy  condition 

L 

material  specimen  model 
thickness 

/  :  wave  interference  condition 

Subscripts : 

n 

number  of  measurements  taken 

cw : 

cooled  wall 

P 

0 

total  pressure 

e  : 

outside  edge  of  boundary  layer 

p 

<1 

static  pressure 

heat  flow 

ex : 

R  : 

jet  tube  exhaust 

radiation 

s 

entropy 

w  : 

heated  wall;  ablation  surface 

T 

temperature 

x  : 

along  the  direction  of  the 

u 

speed  of  gas  flow 

flow  line 

vt 

average  curve  ablation 

*  . 

throat  passage 

speed 

oo  ; 

room  conditions 

X 

length  beginning  from  the 
extension  section  entrance 

1,2: 

before  ablation  and  after 
ablation  or  measurement  points 

1 .  FOREWORD 

When  a  satellite,  a  spacecraft  or  any  such  high  speed  aircraft 
re-enters  the  atmosphere  it  is  going  to  meet  with  extremely  high 
aerodynamic  heating.  Because  of  this,  when  the  various  types  of 
high  speed  aircraft  are  designed,  it  is  necessary  to  carry  out  heat 
protection  calculations.  Experiments  carried  out  on  the  surface  in 
order  to  simulate  material  ablation  confirm  the  theory  of  the  cal¬ 
culations,  give  precise  determinations  of  the  ablation  performance 
of  materials  as  well  as  providing  an  important  foundation  for  beat 
protective  designs  and  effective  measures  to  this  end.  In  the 
re-entry  process  of  spacecraft  and  aircraft,  areas  of  turbulence 
cover  the  greater  part  of  the  surface  of  the  aircraft.  Because  of 
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this,  the  study  of  ablation  under  turbulence  flow  conditions  is 
not  only  absolutely  necessary;  moreover,  the  data  which  is  offered 
by  experiments  in  turbulence  flow  ablation  has  a  very  large  in¬ 
fluence  on  heat  protective  designs. 

At  the  beginning,  people  used  pipe  type  equipment  in  order  to 

carry  out  research  experiments  in  turbulence  flow  ablation  [1], 

Figure  1  shows  the  flow  field  for  a  re-entry  nose  cone  with  a  cone 

shape.  Let  us  assume  that  AB  and  CD  flow  lines  define  the  range 

around  the  outside  of  the  boundary  layers  of  the  object  in  which 

it  is  possible  to  ignore  the  enthalpic  gradient.  These  flow  lines 

on  the  front  surface  of  a  shock  wave  define  the  limits  of  a  ring 

belt  of  which  the  radius  is  R  .  We  recognize  the  fact  that  the 

s 

flow  in  the  ring  belt  on  the  front  circumference  ABCD  of  the  nose 
cone  is  an  equienthalpic  expansion  flow  and  this  flow  will  on  the 
surface  of  the  object  produce  stagnation  point  conditions  at  PQ  and 
Hq.  From  this  condition,  people's  thinking  led  to  the  idea  of 
using  a  jet  tube  with  radius  R,,  to  replace  the  radius  Ro  flow  tube. 
This  would  offer  stagnation  chamber  conditions  at  PQ  and  Hq,  and 
one  could  use  the  jet  tube  flow  to  reproduce  the  equienthalpic 
expansion  flow  around  the  nose  cone  of  the  spacecraft.  Both  inside 
China  and  abroad,  research  workers  employed  pipe  type  equipment  in 
the  performance  of  large  amounts  of  experimental  research  [2,3], 
However,  the  pipe  technique  possesses  certain  inherent  drawbacks. 

The  most  important  of  these  are  the  following.  If  one  uses  this 
type  of  equipment,  observation  is  inconvenient  along  with  the  fact 
that.  In  the  ablation  process,  if  one  increases  the  diameter  of 
the  tube,  the  heat  flow  drops  and  It  is  difficult  to  determine  pre¬ 
cise  values  for  this  kind  of  unstable  effect.  In  order  to  overcome 
these  shortcomings,  workers  in  the  field  carried  out  research  in 
the  area  of  the  experimental  technology  involved  with  turbulent  flow 
flat  plate  ablation. 


2.  BASIC  PRINCIPLES  AND  THE  STRUCTURE  OF  EXPERIMENTAL  EQUIPMENT 
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Concerning  the  high  temperature  air  flow  which  is  shot  out 
from  the  electric  arc  heating  apparatus,  when  it  passes  through 
the  extension  section  it  forms  laminar  flow  boundary  layers  and 
in  the  vicinity  of  the  throat  passage  of  the  jet  tube,  it  turns 
into  a  turbulence  boundary  layer  flow.  At  the  place  where  there 
is  an  exhaust  in  the  wall  of  the  rectangular  supersonic  two-dimen¬ 
sional  jet  tube,  there  is  positioned  a  piece  of  equipment  of  the 
one-flat-plate-model-type  which  is  set  at  a  certain  fixed  angle 
of  attack  to  the  air  flow;  these  two  are  fitted  so  closely  together 
that  there  is  no  gap  between  them  (see  Figure  2).  The  boundary 
layers  which  are  on  this  model  are  natural  extensions  of  the  bound¬ 
ary  layers  on  the  surfaces  of  the  walls  in  the  jet  tube.  From 
these,  on  the  flat  plate  model,  one  obtains  a  fully  developed  tur¬ 
bulent  flow  boundary  layer  and  because  of  the  action  of  the  back 
pressure  gradient  which  is  formed  due  to  the  oblique  shock  wave  of 
the  forward  edge  of  the  flat  plate  model,  one  sees  an  increase  in 
the  simulation  range  of  parameters  in  the  model. 

Ablation  is  an  extremely  complex  process.  It  Is  a  phenomenon 
which  Is  the  result  of  the  combined  actions  of  such  various  factors 
as  thermal  energy  exchange  between  air  flows,  pressure,  shear, 
chemical  reactions  and  so  on.  In  the  ablation  process,  the  related 
parameters  can  be  described  using  the  equations  set  out  below: 


t 

p.  (£*/>)  CO 

-  l(PoJlo)  O,  N,  NO  %3  (2) 

4'~  -  O.iO'Jr.p*,^ /•*)-*>( Re? )-'■ Xu,  -  o  (3) 

r  -  0.82*5^. O  (4) 


1-  or  2-  is  3-  etc. 

The  physical  parameters  p*,  p*  as  well  as  other  similar  para¬ 
meters  are  determined  by  the  equation  below  with  reference  to 
enthalpy  and  on  the  basis  of  calculations  concerned  with  enthalpy, 
that  is 

//*  -  a.::n„  +  o.5A,„  +  n.2M,  (5) 
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The  effective  heat  of  ablation  is  an  important  parameter  to  the 
evaluation  of  the  heat  resistant  performance  of  ablation  mater¬ 
ials.  This  parameter  represents  the  amount  of  heat  which  a  unit 
of  mass  of  an  ablation  material  is  capable  of  absorbing  or  block¬ 
ing  and  its  definition  is  as  follows: 

n.„  -  l <#..(//  -  O./f/o  -  -  E*ri\/pv,  (6 ) 

From  the  various  equations  set  out  above,  the  only  independent 

parameters  there  are  P  or  (p),  H  and  u  .  When  one  is  dealing 

with  a  situation  in  which  the  dimensions  and  configuration  of  a 

re-entry  body  are  determined,  then  the  three  parameters  PQ  or  (p), 

H  and  u  are  only  determined  by  re-entry  speed  (or  M)  and  the 
o  c 

re-entry  altitude.  Because  of  this,  in  these  experiments,  the 

simulation  of  the  three  parameters  PQ  or  (p),  Hq  and  uc ,  is  simply 

based  on  the  reconstruction  of  the  re-entry  ablation  environment. 

To  summarize  what  was  said  above,  in  turbulent  flow  flat  plate 

ablation  experiments,  the  main  parameters  which  must  be  measured 

are  H  ,  P  (or  p),  q  ,  V..  ,  T  and  t  as  well  as  other  related  para- 
o  o  ncw  1  w 

meters . 

The  complete  array  of  experimental  equipment  includes  an  elec¬ 
tric  arc  heating  device,  an  induction  section,  an  extension  sec¬ 
tion  and  rectangular  form  jet  tube,  model  and  measuring  device, 
supports,  water  vapor  and  its  control  system  as  well  as  various 
types  of  testing  systems  and  other  related  gear.  The  various  types 
of  components  involved  were  designed  and  manufactured  on  the  basis 
of  such  principles  as  the  assumption  that  the  power  from  the  elec¬ 
tric  arc  in  the  electric  arc  heating  device,  its  operational  con¬ 
figuration  and  the  oblique  shock  wave  which  is  produced  by  the  main¬ 
tenance  of  the  front  edge  of  the  flat  plate  model  and  the  expansion 
wave  which  is  reflected  back  when  the  oblique  shock  wave  meets  the 
free  boundary  do  not  fall  on  the  model.  In  the  case  of  the  jet 
tube,  a  simple  form  of  analysis  is  employed  and  we  make  use  of  com¬ 
puters  to  carry  out  the  surface  coordinate  calculations.  The  model 
and  the  Jet  tube  involved  have  the  same  width.  The  measuring 


Figure  1.  Spherical-  Figure  2.  Principles  behind  the 

conical  re-entry  nose  experimental  equipment, 

cone  flow  field.  Key:  1 — electric  arc  heating  device; 

Key:  1 — bow-shaped  shock  2 — induction  section;  3 — cooling 

wave  water;  4 — extension  section;  5 — 

jet  boundary;  6 — expansion  wave; 

7 — model  and  support  (or  testing 
device);  8 — oblique  shock  wave 

devices  and  the  ablation  samples  have  heat  absorbing  surfaces  with 
the  same  dimensions.  The  most  important  among  these  are  the  heat 
sink  type  flat  plate  pressure  distribution  measuring  device  (see 
Plate  2(a))  and  the  heat  flow  distribution  flat  plate  measuring 
device  (see  Plate  2(b))  which  are  used  to  make  precise  measure¬ 
ments  of  the  ablation  environment  of  the  sample  models.  The  sup¬ 
ports  are  used  to  hold  the  models  and  the  measuring  devices. 

During  these  experiments,  the  angle  of  deflection  which  was  chosen 
for  the  flow  over  the  models  was  10°. 

3.  PARAMETER  MEASUREMENTS  AND  CALCULATIONS 

In  order  to  obtain  an  ideally  uniform  flow  field,  it  is 
necessary  during  adjustments  to  the  configurations  to  maintain  the 
Jet  tube  exhaust  static  pressure  at  a  value  which  is  basically 
equal  to  the  pressure  of  the  surrounding  environment.  Once  the 
adjustments  have  been  made,  it  is  also  necessary  to  make  use  of 
the  experimental  configurations  set  out  in  Table  1.  In  these  exper 
iments,  we  made  use  of  principles  which  are  essentially  the  same 
as  those  which  were  presented  by  Winovich  in  his  total  pressure 


flow  amount  method  [4],  in  order  to  determine  precise  values  for 
the  total  enthalpy  Hq  in  the  throat  portion  of  the  Jet  tubes  being 
considered.  When  the  amounts  of  flow  are  maintained  precisely 
unchanged  under  conditions  of  hot  and  cold  (by  the  use  of  an 
installed  jet  nozzle  speed  control),  one  makes  use  of  a  cold  air 
standard  to  precisely  determine  the  effective  area  of  the  throat 
passage,  that  is 

G.y/TZ/i96.5P9  (7) 

We  made  use  of  a  turbine  wheel  flow  amount  meter  and  a  heat  sensi¬ 
tive  electrical  resistance  temperature  gauge  in  order  to  make  re¬ 
peated  measurements  of  the  amounts  of  flow  of  the  jet  tube  cooling 
water  and  the  rise  in  temperature  associated  with  this  water.  From 
the  energy  equilibrium  principle,  it  is  possible  to  obtain  the 
total  enthalpy  for  jet  tube  exhausts,  that  is 

-  Ho*  “  aC,Gwh.T/G.  (8  * 


TABLE  1.  Experimental  configuration  parameter  table. 


* 

\3tR& 

f'  \ 

l 

2 

3 

3 

5 

6 

7 

aif/lui*' 

12.  3 

15.1 

10.  3 

10.2 

10.5 

17.1 

17.3 

Ho.. 

7XtlW r 

520 

997 

1500 

2030 

2330 

1336 

1897 

- 

9  «  m  «  * 

/em  m  *  « 

Key:  1 — parameter;  2 — unit;  3 — configuration;  4 — total  gas  flow 
pressure;  5 — kg/cm2;  6 — jet  tube  exhaust  total  enthalpy; 

7 — kcal/kg;  8 — type  of  jet  tube;  9 — conical  surface  jet 
tube;  10 — rectangular  jet  tube 


We  made  use  in  these  experiments  of  heat  sink  type  pressure 
measuring  flat  plates  and  pressure  change  transmitting  devices  to 
measure  the  pressure  distribution  at  five  points  on  the  center 
line  of  the  models  involved;  moreover,  at  the  same  time,  we  took 
measurements  of  the  static  pressure  in  the  Jet  tube  exhaust.  The 
cooling  water  heat  equilibrium  type  heat  measuring  device  (water 
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calorimeter)  measured  the  cold  wall  heat  flow  at  three  points  on 
the  center  line  of  the  flat  plate  model.  In  order  to  make  it 
easier  to  test  produce  a  aero  point  thermal  measuring  device  and 
make  use  of  sweep  measurement  [5]  as  a  substitute  for  fixed  type 
measurement.  We  used  this  approach  to  get  the  heat  flow  distri¬ 
bution  results  which  are  presented  here  (see  Table  2).  During  the 
process  of  ablation  of  a  sample  model  we  made  use  of  an  improved 
photoelectric  colorometric  high  temperature  measuring  device  to 
measure  the  colorometric  temperatures  on  the  ablation  surface  and 
we  recognized  the  fact  that  the  colorometric  temperature  T,***Tm. 
We  also  made  use  of  a  radiation  type  high  temperature  measuring 
device  in  order  to  measure  the  radiated  surface  temperature  of  the 
model  involved,  that  is  T^.  On  the  basis  of  the  Steven-Boltzmann 
principle,  we  can  know  that 

e-(  TjT„y  C9) 


As  far  as  ablation  speed  was  concerned,  we  made  general  use  of 
the  average  time  method  to  determine  it  and  we  recognized  the  fact 
that  ablation  is  of  a  quasi-fixed  nature  and  for  several  models  of 
combustion,  the  ablation  measurements  before  and  after  the  process 
in  terms  of  the  thickness  of  che  models  involved  can  be  obtained 
from  the  equation  below,  that  Is 

a 

Pi  "•  *  2  (*•«  ”  *<«)/'«  (10) 

Concerning  a  few  localized  areas  of  ablation,  if  one  measure  the 
ablation  In  step  and  column  and  other  similar  types  of  ambient 
stall  spots  by  the  use  cf  the  instantaneous  method  of  determination 
one  can  find  the  curve  which  relates  the  amount  of  ablation  A6  to 
the  ablation  time  t;  moreover,  by  solving  for  the  slope  of  the 
curve  at  various  different  instants,  it  is  possible  to  obtain  the 
instantaneous  speed  of  ablation,  that  is 

CL1) 


V„  -  W(A 8),/dt 


Testing  procedures  have  demonstrated  that  the  pressure  p^  at 
the  first  measurement  point  on  the  flat  plate  is  the  static  press¬ 
ure  behind  the  oblique  shock  involved.  We  already  know  that 
S"“P i/f«.  and  if  we  select  the  specific  heat  k  =  1.2,  then  we  get 
[6] 

l)f  +  (*-  1)]  —2(sel -!)}/{*[(«-  1)|  +  (*  +  1)]}  (12) 

When  M1  and  p^  are  precisely  determined,  then  we  can  make  use  of 
the  relationships  presented  below,  that  is 

Mi  -  ujax  -  y/ 2](.U*t*  ~  *.) /«i  (13) 

«i  “  f\Cp\t  ^i)  (1^ ) 

By  repeatedly  consulting  the  thermodynamic  function  table  for  air 
[7],  we  managed  to  directly  obtain  a  simultaneous  solution  for  the 
two  equations  presented  above.  Because  of  this,  it  is  possible  to 
precisely  determine  all  such  parameters  as  a-^ ,  h^,  ug^,  ,  s^  for 
the  point  involved.  We  recognize  the  fact  that  the  gas  flow  over 
the  flat  plate  is  equienthalpic  and  because  of  this  fact,  all  the 
parameters  for  the  various  points  on  the  center  lines  of  the  flat 
plate  models  involved  can  be  solved  for  by  using  the  methods  pre¬ 
sented  above. 


TABLE  2.  Typical  parameter  measurement  results 


Key:  1 — state;  2 — measurement  point;  3 — parameter; 

If — pressure;  5 — water  calorimeter  heat  flow;  6 — zero 
point  heat  flow;  7 — theoretical  heat  flow;  8 —  kilograms 
per  centimeter  square;  9 — calories  per  square  per  second. 
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4.  ANALYSIS  OF  RESULTS  AND  DISCUSSION 

1.  The  distinguishing  characteristics  of  turbulent  flow 
ablation.  After  going  through  experiments  concerning  three  aspects 
of  operation,  we  were  able  to  distinguish  the  various  characteris¬ 
tics  of  turbulence  flow  ablation.  Teflon  (polytetrofluoroehthylene ) 
is  an  excellent  ablation  material.  During  ablation,  the  surface 
temperature  of  it  is  very  low  and  it  reacts  differently  to  laminar 
and  turbulent  flows.  Its  effective  ablation  heat  is  represented 
by  the  expression: 

«...  -4.JW1  (.15) 

Under  five  relatively  high  configurations,  the  ablation  measurements 
which  resulted  for  teflon  are  shown  in  Figure  3  and  these  results 
agree  very  well  with  the  turbulence  flow  ablation  results  found  in 
[6J. 

Besides  this,  by  definition  we  have  the  fact  that  R e, —  V 
and  it  is  shown  by  the  calculations  involving  the  first  measurement 

X  lo7 

point  on  the  flat  plate,  that  Re,  —  •M  x  10'— 4.2^,  and  all  these  values 
are  larger  than  the  critical  turning  point  for  Reynolds  numbers 
5x10 J .  We  made  use  of  formula  (.3)  and  the  Eckert  flat  plate  refer¬ 
ence  enthalpy  method  in  order  to  calculate  the  theoretical  turbu¬ 
lent  flow  heat  flow  distribution  and  the  measured  values  agree  with 
the  corresponding  theoretical  values  within  a  range  of  ±152.  As 
far  as  conical  surface  jet  tubes  are  concerned,  we  made  use  of  the 
equation  below  in  order  to  correct  wave  system  interference  [8], 
that  is 

4tm/4cm*  “  (p'/py*  (16) 

The  results  of  these  calculations  are  set  out  in  Table  2.  The  mea¬ 
sured  values  and  the  theoretical  values  obtained  on  the  basis  of 
corrections  for  wave  system  interference  tend  to  agree  with  each 
other.  Their  numerical  values  are  close  and  the  largest  error 
involved  is  smaller  than  ±172. 


2.  The  influence  of  flow  fields.  Concerning  the  case  in 
which  the  exhaust  flow  field  of  a  conical  surface  jet  tube  is  not 
sufficiently  ideal,  one  can  see  from  the  two  photographs  in  Plate 
2c  which  were  taken  from  orthogonal  directions  that  there  are  both 
expansion  and  compression  wave  systems  in  existence.  However,  if 
one  takes  a  look  at  the  ablation  results,  one  can  see  that  the  wave 
system  interference  has  still  not  created  obvious  troughs  in  the 
sample  model  surfaces.  This  means  that  the  wave  system  interfer¬ 
ence  is  very  weak.  It  can  be  clearly  seen  from  the  measurement 
results  in  Table  2  that  the  surface  jet  tube  flow  fields  involved 
show  a  relatively  large  improvement  and  the  flow  on  the  flat  plate 
closely  approach  the  characteristics  of  equienthalpic  expansion 
waves . 


3.  Ablation  measurement  results.  The  teflon  ablation  measure¬ 
ment  results  are  set  out  in  Table  3.  Plate  2d  is  a  photograph  of 
the  surface  situation  after  teflon  flat  plate  model  ablation  has 
taken  place. 

TABLE  3.  Measurement  results  for  teflon  ablation 


C 

- , - —T - - 

Vi  «*/«> 

l 

2 

3 

6 

5 

0.60 

0.52 

0.51 

0.51 

0.58 

0.68 

0.6-1 

0.59 

0.56 

0.58 

s 

0.78 

0.73 

0.70 

0.68 

0.82 

6 

0.76 

0.71 

0.66 

0.66 

0.63 

7 

0.99 

0.91 

0.89 

0.83 

0.68 

Key:  1 — parameter;  2 — measurement  point;  3 — configuration; 

4 — ablation  speed;  5 — mm/sec 

4,  The  influence  of  a  convex  base  on  ablation.  Due  to  a 
step  effect  in  front  of  a  step  involved,  the  air  flow  stalls  and 
forms  a  stall  shock  wave  (see  Plate  2e).  In  the  stall  spots 
Involved,  heat  flow  and  increases  in  stress  cause  local 


v 
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Figure  3.  A  comparison  of  teflon  ablation  measurement 
results . 

Key:  1 — effective  ablation  beat;  2 — cal/gram;  3 — cal/gram; 

4 — laminar  flow;  5 — turbulent  flow;  6 — total  enthalpy; 

7 — range  of  measurement  values  in  Experiment  I  in 
reference  [4]. 

intensifications  of  ablation.  However,  due  to  the  fact  that  in 
the  ablation  process,  the  steps  are  turned  very  rapidly  into  slopes 
by  the  combustion  (see  Plate  2f).  If  one  adds  to  this  fact  a  case 
in  which  the  Mach  number  of  the  air  flow  involved  is  not  high,  then 
the  effects  of  the  steps  are  not  significant  and  in  front  of  the 
steps  there  is  no  obvious  ablation  trough.  However,  when  one  con¬ 
siders  a  sample  going  through  ablation  measurements  without  steps 
on  the  specimen,  one  can  still  get  the  appearance  of  such  a  trough. 
If  one  is  speaking  in  terms  of  certain  kinds  of  materials,  ablation 
speeds  in  the  area  in  front  of  steps  can  be  one-and-a-balf  times 
higher  than  the  speeds  for  the  same  samples  without  steps.  This 
sort  of  result  is  in  basic  agreement  with  the  results  of  the  calcu¬ 
lations  done  in  [93.  The  results  of  step  model  ablation  have  pro¬ 
vided  precious  experimental  data  with  which  to  further  our  under¬ 
standing  of  the  certain  amount  of  sudden  ambient  ablation  which  is 
unavoidable  with  re-entering  spacecraft  and  aircraft. 

Concerning  the  employing  of  the  technology  of  turbulent  flow 
flat  plate  ablation  experimentation,  there  has  been  research  done 
in  the  analysis  of  various  other  aspects  of  the  subject,  such  as 
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carbon  layers,  heat  decomposition  layer  thicknesses  and  the  nature 
of  these  layers,  the  way  in  which  surface  happens  and  the  pheno¬ 
mena  involved,  insulating  properties,  the  presence  of  mechanical 
erosion  and  the  factors  influencing  it,  the  good  points  and  bad 
points  involved  in  the  working  of  different  types  of  materials  as 
well  as  other  questions  of  a  similar  type.  However,  these  matters 
will  not  be  discussed  here. 

5.  CONCLUSION 

The  turbulent  flow  flat  plate  ablation  experimentation  tech¬ 
nology  has  overcome  the  limitations  of  the  pipe  type  methods.  This 
technology  is  capable  of  quite  adequately  simulating  the  supersonic 
turbulent  flow  ablation  environment  present  on  spacecraft  and  air¬ 
craft  which  are  re-entering  the  atmosphere  and  in  doing  this,  has 
promoted  research  into  the  experimental  technology  related  to  tur¬ 
bulent  flow  ablation.  Such  experimentation  proves  that  this  type 
of  technology  is  an  effective  means  of  running  tests  on  materials 
which  are  heat  resistant  in  terms  of  turbulent  flow  ablation  and 
this  sort  of  technique  expands  the  range  of  experimental  tasks  which 
can  be  undertaken  by  electric  arc  beating  devices. 

What  finally  needs  to  be  pointed  out  is  the  fact  that  the  work 
Involved  in  such  tasks  as  the  measuring  of  enthalpy  by  the  total 
pressure  flow  amount  method,  the  zero  point  measurement  of  heat  flow 
as  well  as  of  surface  temperature,  the  measurement  of  radiant  temper 
ature  and  so  on,  was  completed  by  Wu  Zong-sen,  Li  Ding-wan,  Huang  Yi 
Gu  Shi-jie  and  other  comrades,  and  I  wish  to  express  thanks  for 
this  help. 
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A  PRIMARY  RESEARCH  ON  ABLATION  EXPERIMENT  TECH¬ 
NIQUE  OF  THE  TURBULENT  FLOW  OVER  A 
FLAT  PLATE  SPECIMEN 

If;ui  Yin-da  lVang  Kc-xiang  Yu  Ouang-rotig 

( Chinese  Aerodynamic  Kescnch  and  Development  Center  (CABDC)) 


Abstract 

A  new  method  utlizing  a  flat  plate  (specimen  placed  at  a  small  angle  to  a  super¬ 
sonic  art.  heated  flow  has  been  developed.  This  method  can  surmount  the  disadvantages 
••f  the  usual  pipe  test  method.  Its  observable  recession  rate  is  relatively  uniform,  and 
surface  temperature  and  radiant  heat  flux  can  be  easily  performed. 

Th/s  paper  presents  the  experiment  principles,  the  configuration  of  the  facility, 
measurement  and  theoretical  calculating  method  of  the  parameters.  The  study  was 
accomplished  with  measured  reservoir  pressures  ranging  from  10.2  to  17.3  kg/cm* 

.  ,0,al  eMtU"lpics  trom  520  to  2440  keal/kg.  Under  five  turbulent  conditions  the  abla¬ 

tion  t»st  results  of  teflou  models  were  given  in  this  paper. 
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A  NEW  METHOD  OF  COOLING  TURBINE  VANES 


Gu  Wei-sao,  Zhang  Yu-ming  and  Xu  Hong-kun 

(Institute  of  Engineering  Tbermopbysics ,  Academia  Sinica) 

As  far  as  current,  state-of-the-art  turbine  blades  in  avia¬ 
tion  engines  are  concerned,  the  generally  employed  cooling  method 
is  a  complex  one  involving  blast,  convection  and  gas  films;  in 
these  same  advanced  blades,  the  complexity  of  the  construction  is 
getting  more  and  more  extreme  and  the  amounts  of  cooling  gas  in¬ 
volved  is  getting  larger  and  larger.  In  order  to  raise  still  fur¬ 
ther  the  intake  temperatures  of  turbines,  it  is  necessary  for  one 
to  look  for  a  new  way  of  handling  cooling.  This  article  suggests 
the  employment  of  an  internal  cooling  process  within  strengthened 
blades  along  with  the  employment  of  additives  to  lessen  the  trans¬ 
fer  of  heat  to  the  blades  from  the  fuel  gases.  The  combination  of 
these  two  measures  would  improve  cooling.  In  the  article  below, 
the  authors  deal  respectively  with  the  bases  for  these  measures  and 
the  temperature  calculations  involved  with  regard  to  the  first 
stage  guide  vanes  of  the  MK202  engine. 

1.  COARSENESS  STRENGTHENED  HEAT  EXCHANGE  IN  THE  INNER  PASSAGEWAYS 
OF  BLADES 

In  engineering,  the  use  of  coarse  surfaces  to  strengthen  heat 
exchange  is  already  widely  in  use  and  there  are  some  advanced  avia¬ 
tion  engines  which  have  blade  cooling  systems  which  also  employ 
coarse  surface  inner  cooling  structures.  The  basic  principle  in¬ 
volved  in  coarseness  strengthened  heat  exchange  is  that  the  coarse¬ 
ness  interferes  with  the  air  flow  and  disrupts  lower  layer  viscosity. 
This,  as  a  result,  lowers  the  resistance  to  beat  exchange  between 
the  wall  surfaces  and  the  cooling  gases.  Generally,  the  level  of 


This  article  was  read  to  the  3rd  Annual  All-China  Engineering 
Thermophysics  Technology  Conference  in  Qui  Zhou  in  April  of  1980. 


the  strengthening  of  the  heat  exchange 
is  related  to  the  configuration  of  the 
coarseness,  the  dimensions  of  it,  the 
geometry  of  the  flow  path,  the  value 
of  the  flow  body  Pr  and  the  flow  Re 
value.  What  is  shown  in  Figure  1  are 
the  curves  defined  by  the  relationship 
between  the  heat  transfer  strength 
ratio  Nu/Nu  and  the  corresponding 
heights  h/d  nodal  height  ratios  s/h. 
One  only  needs  to  select  appropriate 


<  h 


Re  =  10*,  Pr  =  0.71 


Figure  1.  Heat  transfer 
strengthening  ratios  for 
coarsened  cylinders 


parameters  of  coarseness  and  it  then  becomes  possible  to  effect¬ 
ively  raise  the  heat  transfer  coefficient  for  the  inner  cavities 
of  the  blades  concerned. 


The  cooling  structures  inside  of  blades  usually  take  the  form 
of  unusually  shaped  passageways  and  in  the  acute  angle  sections  of 
these  passageways,  it  is  possible  to  put  the  air  flow  into  a  lam¬ 
inar  flow  configuration.  In  such  circumstances,  the  coarsened 
support  surfaces  have  no  strengthening  effect.  Because  of  this, 
it  is  possible  to  make  use  of  horizontal  supports  in  the  vicinity 
of  even  a~eas ,  depending  on  flow  area  interference  with  the  neigh¬ 
boring  supports.  At  the  same  time,  in  the  end  sections  of  the  angled 
areas,  if  one  adds  a  concave  trough  in  order  to  increase  the  degree 
of  turbulence  of  the  flow  in  the  angled  area,  then  one  can  guarantee 
the  required  strength  of  heat  transfer. 

2.  THE  EMPLOYMENT  OF  ADDITIVES  TO  LOWER  THE  TRANSFER  OF  HEAT  FROM 
THE  FUEL  TO  THE  BLADES 


After  the  ;ombustion  of  a  fuel  that  contains  small  amounts  of 
(SiOCgHS)^,  the  SiO^  which  is  produced  in  fuel  gases  whose  temper¬ 
atures  are  lower  than  1700K  will  exist  in  the  form  of  suspended 
particles  with  diameters  of  0.1-1  microns.  The  amount  of  these 
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particles  which  passes  through  a  certain  integral  volume  surface 
of  flow  can  be  represented  as 

J  “*  Jp  +  J »■  4-  J,  •=  — /),  grad  n  +  Wn  +  V,/)  ) 

In  this  equation  J^,  and  J^,  respectively  represent  the  amount 
of  molecular  diffusion  of  particles,  the  amount  of  convection 
diffusion  of  particles  and  the  amount  of  particles  which  are  caused 
to  move  from  a  relatively  higher  temperature  to  a  relatively  lower 
temperature  due  to  the  effects  of  thermal  drag  in  uneven  tempera¬ 
ture  fields.  is  the  coefficient  of  diffusion  of  these  particles 

in  air.  n  is  the  concentration  of  particles.  VT  is  the  speed  of 
movement  of  the  particles  relative  to  the  gas  involved  as  induced 
by  thermal  drag. 

The  most  important  influence  in  the  core  of  a  turbulent  flow 
is  the  amount  of  convection  diffusion  J^.  When  particles  follow 
the  same  air  flows,  they  quickly  diffuse.  However,  in  viscous  low 
layers  Jw  is  very  close  to  zero.  To  put  it  the  other  way  around, 

'  JT  In  the  core  of  a  turbulent  flow  can  be  Ignored;  however,  in  a 

viscous  low  layer.  It  becomes  the  deciding  factor  which  causes  the 
particles  to  precipitate  out  onto  a  cooled  wall  surface.  According 
to  the  Epstein  formula  [3] 

Vt  *“  pT )  grad  T  (2) 

In  viscous  lower  layers  i ■“  -*i(.3P/3y)y<«  “  “(T*  *"  ,  From 

this  one  can  obtair.(a7’/3y),<j  *“  “C7”*  —  and  after  one  substitutes 

this  Into  Equation  (2),  because  of  the  fact  that,  as  far  as  the  fuel 
gases  are  concerned,  0/i)><4  (ft  i)g,(pT),<t  *^PtTt>  ,  It  is  finally 

possible  to  obtain  the  viscous  lower  layer  precipitation  speed  for 
SiOg  particles,  which  is 

Vjw - /?(aPf/pC,)f(  1  -  TJT, )  (  3  ) 

The  negative  sign  represents  the  direction  of  VTW  which  is  opposite 
to  y.  It  can  be  seen  from  this  that  under  similar  conditions  of 

I 

temperature  and  pressure,  the  smaller  the  ratio  value  is  between 
the  wall  temperature  T  and  the  gas  flow  temperature  T  ,  or  the 

W  g 
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Figure  2.  Heat  Flow  Drop 
in'  Throat  Section  of  Jet  Tube 

1  -  Si02  content  ( % ) 

larger  the  coefficient  of  heat  transfer  is,  the  greater  will  be 
the  speed  of  precipitation  of  Si02  particles  on  the  cool  wall 
surfaces . 

The  use  of  additives  to  lower  the  jet  tube  heat  flow  in  -rock¬ 
et  engines  has  already  obtained  practical  results.  For  example, 
the  U.  S.  Agena  lifting  rocket  uses  a  1%  addition  of  a  silicone 
compound  which  causes  on  the  walls  of  its  jet  tubes  the  formation 
of  an  active  covering  layer  of  Si02  and  this  successfully  lowers 
the  flow  of  heat  [4],  Figure  2  displays  results  of  the  measurement 
of  heat  flow  in  the  throat  sections  of  jet  tubes  as  these  measure¬ 
ments  were  done  in  [5].  Looking  from  the  point  of  view  of  research 
materials  available  from  inside  China  and  abroad,  the  addition  of 

small  amounts  of  additives  to  fuel  can  cause  jet  tube  heat  flow  to 

SiO, 

be  lowered  in  the  neighborhood  of  30-40J.  The^ precipitation  layer 
thickness  in  various  sections  of  the  jet  tubes  involved  runs 
approximately  from  0.03-0. 17mm.  Obviously,  the  larger  the  coeffi¬ 
cient  of  heat  exchange  is,  the  larger  will  be  the  effect  which  is 
caused  by  the  heat  resistance  of  the  Si02  precipitation  layer  on 
the  lowering  of  the  jet  tube  heat  flow. 


In  order  to  do  research  on  the  effects  which  fuels  containing 
carbon  have  on  the  cooling  performance  of  high  temperature  turbine 
blades,  [6]  added  a  small  amount  of  ethylene  silicate  to  kerosene 
fuel  and  after  combustion  of  this  mixture,  obtained  temperatures 
of  1300-1900K.  During  experiments  in  which  fuel  containing  Si02 
passed  through  blade  cascades  of  cooling  vane  blades,  there  was  a 
precipitation  layer  of  Si02  particles  formed  on  the  wall  surfaces 
of  the  blades  and  this  layer  had  a  uniform  thickness  of  0.15mm. 
Because  of  this,  experimentation  and  theory  both  demonstrate  that 
gas  flows  which  have  small  amounts  of  Si02  particles  suspended  in 
them  are  capable  of  forming  a  thin  active  covering  layer  on  the 
cooling  surfaces.  This  increases  heat  resistance  and  because  of 
this,  it  is  possible  to  reduce  the  transfer  of  heat  from  the  fuel 
gas  to  the  blades. 

3.  CALCULATIONS  FOR  THE  IMPROVEMENT  OF  THE  DESIGN  OF  FIRST  STAGE 
GUIDE  VANES  IN  MK202  ENGINES 

On  the  basis  of  the  design  materials  available  concerning  the 
first  stage  guide  vanes  of  MK202  Spey  engines  [7],  if  one  takes  the 
turbine  wheel  intake  temperature  and  raises  it  to  1600K,  then  the 
corresponding  peak  fuel  temperature  value  is  1920K;  the  temperature 
and  the  amounts  of  flow  of  the  cooling  gases  remain  constant  at 
78OK  and  2.9$.  In  order  to  strengthen  the  interior  cooling  process 
in  turbine  blades,  the  flow  path  is  changed  to  the  structure  which 
is  shown  in  Figure  3.  The  cooling  air  enters  from  the  base  section 
of  the  forward  cavity,  then  passes  through  the  forward  cavity,  the 
middle  cavity  and  the  after  cavity  by  way  of  the  finned  passages 
through  them  and  finally,  passes  through  a  row  of  0.6mm  holes  in 
the  tail  section  of  the  blades  concerned  and  it  is  expelled  through 
the  dish  of  the  blades  involved.  On  the  interior  wall  surfaces  of 
the  blade  backs  and  blade  dishes  in  the  cooling  passageways,  there 
are  cast  horizontal  fins  with  coarse  surfaces.  The  dimensions  of 
these  surfaces  uniformly  correspond  to  the  conditions: 
h/d  ■■  0.02,  s/h  “•  10  ,  and  the  coefficient  of  heat  exchange  involved  is 
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Figure  3.  The  internal  cooling  structure  of  guide  vanes  and  the 
temperature  distribution  (°C)  as  calculated  on  the  basis  of  a 
fuel  gas  temperature  of  1920K 

calculated  on  the  basis  of  the  Webb  heat  transfer  formula  coarse 
round  tubes.  Concerning  the  temperature  rise  in  the  cooling  gases 
in  the  interior  of  the  various  passages,  it  is  possible  to  assume 
that  the  interior  wall  temperatures  along  constant  heights  do  not 
vary  and  that  they  are  equal  to  the  average  value  Tw2i>  From 
these  assumptions,  one  can  deduce  the  calculation  formula  that 
follows,  that  is 

-  ( TwU  -  T'clX  1  -  *-)  (4  1 

In  this  equation  m  —  aeFjGcc,„  T'„  is  the  cooling  gas  intake  temp¬ 
erature  for  the  passage  under  consideration.  The  external  heat 
transfer  coefficient  for  the  blades  involved  are  selected  on  the 
basis  of  the  sections  presented  in  the  material  found  in  [7J.  The 
forward  station  is  8710  and,  as  a  consequence,  the  blade  backs  and 
blade  dishes  involved  have  their  sections  respectively  selected  as 
the  average  values  shown  in  Figure  3.  The  calculations  demonstrate 
that  the  highest  temperatures  in  the  tail  sections  of  the  blades 
start  where  the  exhaust  holes  on  the  blade  backs  are  located. 

Because  of  the  fact  that  this  place  is  the  point  at  which  the  exter¬ 
nal  heat  transfer  coefficient  is  at  its  highest  in  the  section 

2 

being  considered,  this  coefficient  is  equal  to  4725  kcal/m  hr  °C. 

In  this  calculation,  in  the  area  left  over,  take  the  thickness  of 
the  Si02  precipitation  layer  to  be  0.05mm  and  the  coefficient  of 
heat  conductivity  to  be  0.53  kcal/m  hr  °C  on  the  basis  of  materials 
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being  used  routinely  outside  of  China. 

Besides  this,  in  order  to  improve  the  cooling  of  the  forward 
edges  of  the  blades  concerned,  make  the  radius  of  the  interior 
arc  of  the  forward  edges  equal  to  the  external  arc.  This  will 
make  it  easier  to  cause  the  interior  and  exterior  heat  trans¬ 
mission  surfaces  equal  at  the  forward  stations.  It  will  also  help 
to  increase  the  forward  edge  thickness  from  0.76  to  1.5mm  and  to 
cast  concave  troughs  with  coarse  surfaces  on  the  inside  walls  of 
the  leading  edges. 

The  blade  temperature  distribution  which  is  calculated  on  the 
basis  of  these  conditions  is  as  shown  in  Figure  3*  With  the 
exception  of  the  tail  sections,  the  temperatures  of  all  the  various 
sections  are  lower  than  the  permissible  operating  temperatures  for 
the  metal  in  the  blades.  If  one  takes  the  amount  of  the  cooling 
gas  and  increases  it  to  (this  is  still  far  lower  than  the  amount 
of  cooling  gas  which  is  required  by  the  standard  cooling  methods, 
that  is,  6-7?)  then  it  is  possible  to  cause  the  maximum  tempera¬ 
tures  for  the  tail  section  to  fall  to  108l°C. 

Finally,  it  still  needs  to  be  pointed  out  that  the  particles 
of  SiC>2  in  the  fuel  gases  will  also  form  precipitation  on  the  sur¬ 
faces  of  the  moving  blades.  This  sort  of  phenomenon  is  capable  of 
greatly  reducing  the  amount  of  cooling  gas  required  by  moving 
blades.  The  rules  which  govern  the  precipitation  of  these  SiC^ 
particles  onto  the  surface  rotor  blades  out  to  be  given  specialized 
study. 
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RESEARCH  OH  THE  RECI/^JLATING  FLOW  RATE  IN  THE  WAKE  OF  TWO- 
DIMENSIONAL  FLAME  STABILIZERS  (UNDER  COLD  CONDITIONS)* 

Zhang  Qing-fan 

(Nanjing  Aeronautical  Institute) 

THEORETICAL  ANALYSIS 

The  recirculating  flow  of  blunt  body  flame  stabilizers  is 
an  important  factor  influencing  flame  stability  characteristics. 

On  the  basis  of  jet  flow  theory  it  is  possible  to  deduce  an  equa¬ 
tion  representing  the  recirculating  flow.  In  Figure  1,  due  to  the 
induction  effects  of  the  jets  involved,  the  amounts  of  flow  grad¬ 
ually  increase.  Let  us  assume  that  the  properties  and  temperatures 
of  the  jets  and  ambient  gases  are  the  same.  In  such  a  case,  if 
we  do  not  figure  in  compressibility,  then  the  total  quantity  of 
flow  in  jets  on  any  plane  a-a  is 

Mm  —  Powoyi  +  po  (  «^y  Cl ) 

Jft 

In  this  equation,  pQ  is  the  jet  density,  u  Is  the  speed 
inside  the  jet  boundary  layer.  On  the  basis  of  these  relationships, 
the  amount  of  induced  jet  flow  within  the  range  x  is 

—  Po“o*o  C2  ) 

In  the  limited  jets  which  are  shown  in  Figure  2,  the  jet  bound¬ 
aries  and  the  wall  surfaces  intersect  at  F  and  the  total  quantity 
of  Jet  flow  on  the  plane  Xp  is 

Mmr  —  prf*»ytT  +  p*  \  "Jy  (3 ) 

In  this  relationship,  y2p  is  half  the  height  of  the  passage  involved. 

I - 
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Thermophysics  Technology  Conference  in  Gui  Zhou  in  April  of  1980. 
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On  the  basis  of  this,  the  largest  amount  of  induced  jet  flow  in 
the  two-dimensional  passage  involved  is 

Mm «  “  M ,T  —  ) 

The  gases  of  the  jet  flow  which  is  induced  by  the  limited 
jets  involved  here  form  in  the  tube  exhaust  a  secondary  flow, 
the  speed  of  which  is  us .  Experiments  demonstrate  that  when  the 
amount  of  flow  involved  in  this  secondary  flow  is  smaller  than 
the  amount  of  induced  flow  in  Equation  (^),  one  will  see  the 
initial  formation  of  a  reflux  flow  first  in  the  vicinity  of  point 
F.  The  smaller  the  amount  of  flow  involved  in  the  secondary  flow 
is,  the  larger  the  area  of  the  reflux  flow  will  be.  When  the 
secondary  flow  is  equal  to  zero,  the  amount  of  induced  jet  flow 
is  equal  to  the  amount  of  reflux  flow.  In  such  a  circumstance, 
the  dimensions  of  the  area  of  reflux  flow  are  at  their  largest, 
filling  the  entire  exterior  part  of  the  passage.  This  Is  equiva¬ 
lent  to  the  exterior  part  of  the  passage  being  blocked  (Figure  3). 
If  one  does  not  figure  in  the  friction  between  the  gas  flow  and 
the  wall  surfaces  involved,  then  the  boundary  conditions  for  Figure 
3  and  Figure  4  are  the  same.  Because  of  this,  one  can  replace 
Figure  3  and  take  the  height  of  the  passage  B  as  half  the  trough 
width  of  the  stabilizer.  It  can  be  seen  that  the  amount  of  flow 
from  the  stabilizer  is  related  to  the  calculated  values  of 
Equation  (4). 

If  the  flow  characteristics  of  limited  jets  and  free  jets  are 
the  same,  then  it  is  possible  to  directly  employ  the  boundary  lines 
presented  in  the  first  section  of  the  related  equations  in  [2], 
that  is 


y,  “  A#  — "  0.4 1 W 

C 51 

The 

external  boundary  line 

y,  At  +  0.584A 

(6) 

for 

the  initial  stage 

The 

Jet  boundary  layer  width 

*  —  y»  —  y\  —  cx 

(7) 
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:'!--ure  1.  A  simplified  diagram 
cf  a  free  planer  jet  flow 


Figure  2.  A  simplified  diagram 
of  limited  planar  jet  flow 


Figure  3.  Planar  Limited 
Jet  When  ug  =  0 


Figure  4 .  Two  Dimensional 
Flame  Stabilizer  Wake  Flow 


In  this  equation,  c  is  a 
coefficient.  The  speed  distribu¬ 
tion  for  the  interior  of  the 
boundary  layers  involved  is 

(«o  —  «)/«0  —  ( 1  —  ijl  5)J  (8 ) 

In  this  equation,  n  is  a  non- 
dimensional  coordinate  and  is 
also  defined  by  the  equation 

n  —  (yt  —  y)/b  (9) 


Figure  5.  The  relationship 
between  the  coefficient  K  and 
the  geometrical  parameters  of 
stabilizer  devices. 

Key:  1 — apex  angle 


One  can  take  the  Equations  ( 5 )  —  ( 9 )■  and  substitute  them  into 
Equation  (3).  Then  one  can  solve  for  the  maximum  amount  of  Induced 
Jet  flow  from  Equation  (4),  that  Is 
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Mt/pohofio  ■*  0.23 B l Mq 


Cio) 


Equation  (10)  represents  the  induced  flow  energy  from  a 
limited  jet  when  the  secondary  flow  is  sufficiently  large  (no 
reflux  flow).  However,  when  the  secondary  flow  is  zero,  the  flow 
conditions  at  the  boundaries  of  the  jet  flow  begin  to  change,  and 
this  is  particularly  true  of  the  blockage  ratio  of  the  stabilizer 
apparatus,  that  is  e— B/(S  +  A0)  .  Therefore,  the  actual  amount 

of  reflux  flow  is  M  ,  or 

Jlf  r/p#A0u0  —  Ke/(  1  —  *)  ^ 

This  is  nothing  else  but  the  equation  representing  the  amount 
of  reflux  flow  as  that  amount  can  be  deduced  from  the  theory  of 
jets.  The  coefficient  K  in  this  equation  can  be  fixed  to  a  precise 
value  by  experimentation. 

METHODS  AND  RESULTS  OF  EXPERIMENTATION 

Concerning  the  measurement  of  the  parameters  of  the  reflux 
flow  area  in  a  cold  condition,  the  dimensions  of  the  transparent 
two-dimensional  passage  are  100x150mm  .  The  length  of  the  stabi¬ 
lizer  device  is  100mm.  The  average  flow  speed  in  front  of  the 
stabilizer  device  is  3^m/s.  By  the  use  of  normal  total  pressure 
values  and  round  dish  investigation  techniques  we  can  respectively 
determine  the  reflux  flow  area  total  static  pressure  and  the  reflux 
flow  introduced  in  the  lower  sections  to  be  the  largest  amounts  of 
reflux  flow  in  the  center  cross  sections.  The  experimental  results 
for  three  types  of  blockage  ratios  and  three  types  of  apex  angles 
for  each  of  nine  types  of  stabilizer  devices  are  presented  in 
Figure  5.  From  this  one  can  obtain  the  K  value  for  Equation  (11), 
which  is 

K-  0.155  +  0.08(0/30°  -  1)  (_12) 

In  this  equation  0  is  the  apex  angle  of  the  stabilizer  device 
(.degrees).  On  the  basis  of  this,  according  to  the  theory  of  jets, 
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one  can  deduce  the  amount  of  reflux  flow  from  a  two-dimensional 
stabilizer  device  in  a  cold  condition,  that  is 

M,/M.  -  [0.155  +  0.08(6/30° 

-l)]e/(l-e)  (13) 

In  this  equation,  M  is  the  total  quantity  of  flow  over  the 
stabilizer  device.  The  appropriate  range  for  use  with  Equation 
(13)  is  this:  stabilizer  device  blockage  ratio  0.2-0.14,  and  stab¬ 
ilizer  device  apex  angle  30-60°.  The  Reynolds  number  defined  by 
the  width  of  the  trough  is  (0 . 7-1 • ^ )xl0^ .  This  article  does  not 
contain  any  research  into  the  influence  of  the  degree  of  turbulence 
in  a  flow  and  because  of  this,  strictly  speaking,  the  results 
obtained  are  only  suitable  for  use  under  experimental  conditions 
in  which  there  are  specially  determined  Reynolds  numbers. 
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A  NEW  METHOD  OF  COOLING  TURBINE  VANES 

On  Wei-zno,  Zhang  Yu-niing,  Xu  Hong-kun 
( hiafitute  of  Kniiincerintf  Thrrmophysics,  Acadimin  Sintra) 

Abstract 

In  this  paper  the  methods  of  augmentation  of  convective  heat  transfer  in  thff 
cooling  channels  of  the  vane  with  rough  surface  and  formation  of  the  dynamic  coating 
on  the  airfoil  surface  in  virtue  of  silicone  additives  in  the  fuel  were  suggested.  Ry 
means  of  these  could  be  raised  the  turbine  inlet  pas  temperature  to  1000  K  when  coolant 
to  gas  flow  ratio  is  equal  to  O.IKI- 41.04. 
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